1. Introduction 


The need to solve optimization and control problems arises in many settings. 
Although in some cases these problems may be easily solved, either analytically or 
computationally, in many other cases substantial difficulties are encountered. For 
example, candidate optimal states and controls may belong to infinite dimensional 
function spaces and one may have to minimize a nonlinear functional of the state and 
control variables subject to nonlinear constraints that take the form of a system of 
partial differential equations whose solutions are in general not unique. In this paper, 
our goal is to construct and analyze a framework for the approximate solution of 
many such problems. The setting for our framework is a class of nonlinear control or 
optimization problems which is general enough to be of use in numerous applications. 
The major steps in the development and analysis of our framework are as follows: 

• define an abstract class of nonlinear control or optimization 
problems; 

• show that, under certain assumptions, optimal solutions exist; 

• show that, under certain additional assumptions, Lagrange mul- 
tipliers exist that may be used to enforce the constraints; 

• use the Lagrange multiplier technique to derive an optimality 
system from which optimal states and controls may be deduced; 

• define algorithms for the approximation, in finite dimensional 
spaces, of optimal states and controls; and 

• derive estimates for the error in the approximations to the op- 
timal states and controls. 

Two of the key ingredients used to carry out the above plan are a theory given in 
[21] for showing the existence of Lagrange multipliers and a theory first developed in 
[6] for the approximation of a class of nonlinear problems. In both of these theories, 
certain properties of compact operators on Banach spaces play a central role. We 
point out that the nonuniqueness of solutions of the nonlinear constraint equations 
deems it appropriate to employ Lagrange multiplier principles. 

After having developed and analyzed the abstract framework, we will apply it 
to some specific, concrete problems. In each case, we use the abstract framework to 
analyze the concrete problems by merely showing that the latter fit into the former. 
The particular applications we consider are: 

• control problems in structural mechanics having geometric non- 
linearities that are governed by the von Karman equations; 

• control problems in superconductivity that are governed by the 
Ginzburg-Landau equations; and 

• control problems for incompressible, viscous flows that are gov- 
erned by the Navier-Stokes equations. 

In considering these applications, we will purposely choose different types of controls 
in order to illustrate how these can be accounted for within the abstract framework. 
In all three cases, approximation will be effected through the use of finite element 
methods. 
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2. The abstract problem and its analysis 


In this section we define and analyze an abstract class of constrained nonlinear 
control problems; an outline of the definitions and results of this section is as follows. 

• In §2.1, the abstract class of constrained control problems that 
we consider is defined. 

• In §2.2, a list of assumptions about the class of abstract prob- 
lems is given. 

• In Theorem 2.1 of §2.3, some of the assumptions listed in §2.2 
are used to show that optimal solutions of the abstract problem 
exist. 

• In §2.4, some additional assumptions of §2.2 are used to show 
that Lagrange multipliers exist that may be used to enforce the 
constraint; also, first-order necessary conditions are given. 

• In Theorems 2.5 and 2.6 of §2.4, the first-order necessary condi- 
tions for determining optimal states and controls are simplified 
under additional assumptions about the control set. 

• In §2.5, the optimality system from which optimal controls and 
states can be determined is made more amenable to approxima- 
tion by simplifying the dependence of the objective functional 
on the control. 

2.1. The abstract setting 

We begin with the definition of the abstract class of nonlinear control or opti- 
mization problems that we study. 

We introduce the spaces and control set as follows. Let G, X , and Y be reflexive 
Banach spaces whose norms are denoted by || • ||g, || • ||x, and || • ||y, respectively. Dual 
spaces will be denoted by (•)*. The duality pairing between X and X* is denoted 
by (• one similarly defines (•, -) Y and The subscripts are often omitted 

whenever there is no chance for confusion. Let 0, the control set, be a closed convex 
subset of G. Let Z be a subspace of Y with a compact imbedding. Note that the 
compactness of the imbedding Z C Y will play an important role. 

We assume that the functional to be minimized takes the form 

(2.1) J(v,z) = XX(v) + X£(z) M{v,z) x 0, 

where T is a functional on X , £ a functional on 0, and A is a given parameter which 
is assumed to belong to a compact interval A C R + . 

The constraint equation M(v,z) = 0 relating the state variable v and the control 
variable z is defined as follow. Let A be a differentiable mapping from X to Y, K a 
continuous linear operator from 0 to Y, and T a continuous linear operator from Y 
to X . For any A £ A, we define the mapping M from X x 0 to X by 

(2.2) M(v, z) = v + XTN(v) + XTK(z) V (v, z) e AT x 0 . 

With these definitions we now consider the constrained minimization problem: 

(2.3) min J(v,z) subject to M(v,z) = 0. 

(v ,z)£X X © 
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In (2.3), we seek a global minimizer with respect to the set {(v,z) £ X x 0 : 
M(v,z) = 0}. Although, under suitable hypotheses, we will show that the problem 
(2.3) has a solution, in practice, one can only characterize local minima, i.e. , points 
(u, g) £ X x 0 such that for some e > 0 

J(u, g ) < J(v, z) V (v, z) £ X x 0 such that 

M(v, z) = 0 and ||m — v\\x < e . 

Thus, when we consider algorithms for locating constrained minima of J , we must be 
content to find local minima in the sense of (2.4). 

After showing that optimal solutions exist and that one is justified in using the 
Lagrange multiplier rule, we will introduce some simplifications in order to render the 
abstract problem (2.3), or (2.4), more amenable to approximation. The first is to only 
consider the control set 0 = G. The second is to only consider Frechet differentiable 
functionals £(■) such that the Frechet derivative £'(g) = E~ 1 g, where E is an invertible 
linear operator from G* to G. 

2.2. Hypotheses concerning the abstract problem 

The first set of hypotheses will be invoked to prove the existence of optimal 
solutions. It is given by: 

(HI) inf^gx T(v) > — oo ; 

(H2) there exist constants ex, j3 > 0 such that £(z) > a||z||^ V z £ 0 ; 

(H3) there exists a (v, z) £ X x 0 satisfying M(v, z) = 0 ; 

(H4) if u in X and g^E g in G where {(«("), g*-"- 1 )} C X x 0, then 

N(u( n ')) N(u) m Y and K(g^) K(g) m Y ; 

(H5) J(-, •) is weakly lower semicontinuous on X x 0 ; and 

(H6) i/{(«( n ),j("))}cXx0is such that {jF(t/( n ))} is a bounded set in IR and 

M(u( n \g(E) = 0, then {«(")} is a bounded set m X. 

The second set of assumptions will be used to justify the use of the Lagrange 
multiplier rule and to derive an optimality system from which optimal states and 
controls may be determined. The second set is given by: 

(H7) for each z £ 0, v J(v, z) and v M(v,z) are Frechet differentiable ; 

(H 8 ) z £(z) is convex, i.e., 

£(jz l + (1 - j) z 2 ) < 7 £( z i) + (1 - 7 )£( z 2 ) V zi, z 2 £ 0 , V 7 £ [0, 1] ; 
and 

(H9) for v £ X, N'(v) maps X into Z. 

In (H9), N' denotes the Frechet derivative of N. 

A simplified optimality system may be obtained if one invokes the additional 
assumption: 

(H10) 0 = G, and the mapping z £( z ) is Frechet differentiable on G. 

Hypotheses (H7)-(H10) allow us to obtain a simplified optimality system for al- 
most all values of the parameter A £ A. In many cases, it is possible to show that 
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the same optimality system holds for all values of A. The following two additional 
assumptions which will only be invoked in case (1/A) is an eigenvalue of —TN'(u) 
each provides a setting in which this last result is valid: 

(Hll) if v* G X* satisfies (i + A [A , (m)]*T*)i>* = 0 and K*T*v* = 0, then 
v* = 0 ; or 

(Hll)' the mapping (v, z) h- > v + A TN'(u)v +A TKz is onto from X x G to Y . 

In order to make the optimality system more amenable to approximation and 
computation, we will invoke the following additional assumption: 

(H12) Z'id) = E _1 9 ; where E is an invertible linear operator from G* to G and 
g is an optimal control for the constrained minimization problem (2.4). 

2.3. Existence of an optimal solution 

We first use assumptions (H1)-(H6) to establish that optimal solutions exist. 

THEOREM 2.1. Assume that the functional J and mapping M defined by (2.1) and 
(2.2), respectively, satisfy the hypotheses (H1)-(H6). Then, there exists a solution to 
the minimization problem (2.3). 

Proof: Assumption (H3) simply asserts that there is at least one element of X x 0 that 
satisfies the constraint. Thus, we may choose a minimizing sequence {(«("), </ n ))} C 
X x 0 such that 

lim J{vf n \ </")) = inf J(v,z) 

00 (v,z) eXx& 

and 

M{u {n \g {n y) = 0 . 

By (HI) and (H2), the boundedness of </ n ))} implies the boundedness of 

the sequences {||<7 ( - n ' ) ||G} and {jF(«( n ))}. Then, by (H6), we deduce that {Hw^/lx} is 
bounded. Thus, we may extract a subsequence {(«("), <j/ n i)} such that i/( n ) u in X 

and </") g in G. Since 0 is closed and convex, we have g G 0. Of course, u G X . 
We next show that (u,g) satisfies the constraint equation. Using (H4), we have that 

lim {TN(u^), f) = lim {N(u^),T*f} = {N(u),T*f} = {TN(u), f) V / G A* 

n— >00 n— >00 

and 

lim (TK(g^)J) = lim (K(g^),T* f) = (K(g),T*f) = (TK(g)J) V / G A* 

n— >00 n— >00 

so that 


0 = lim (M(«("), </(")), /) = (« + A TN(u) + A TK(g), f) V / G A* , 

n— >00 

i.e. , M(u,g) = 0. Finally, we use (H5), the weak lower semicontinuity of J(-, •), to 
conclude that (u, g) is indeed a minimizer in A x 0 satisfying the constraint M(u,g) = 

0 . □ 

2.4. Existence of Lagrange multipliers 

We now wish to use the additional assumptions (H7)-(H9) to show that the La- 
grange multiplier rule may be used to turn the constrained minimization problem 
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(2.3) into an unconstrained one. Actually, the Lagrange multiplier rule will only en- 
able us to find local minima in the sense of (2.4). We first quote the following abstract 
Lagrange multiplier rule whose proof can be found in [21]. 

THEOREM 2.2. Let X\ and X 2 be two Banach spaces and 0 an arbitrary set. Suppose 
J is a functional on X\ x 0 and M a mapping from X\ x 0 to X 2 . Assume that 
( u,g ) E X\ x 0 is a solution to the following constrained minimization problem: 

M(u, g) = 0 and there exists an e > 0 such that J(u,g) < J(v, z) 
for all ( v , 2 ) G X\ x 0 such that ||m — v\\x 1 < e and M(v, z) = 0 . 

Let U be an open neighborhood of u in X\. Assume further that the following conditions 
are satisfied: 

(2.6) for each z G 0, v 1 — > J{v, z) and v 1 — > M(v,z) are Frechet- differentiable at 
v = u ; 

(2.7) for any v G U, 21,22 E 0, and 7 E [0,1], there exists a z 1 = z 1 {v, 21 , 22 ) 
such that 

M(v, z 7 ) = jM(v, z 1 ) + (1 - 7 )M(v, z 2 ) 

and 

J(v, z 7 ) < 7 J(v, z 1 ) + (1 - j)J (v, z 2 ) ; 

and 

(2.8) Range(M [1 ( «,(/)) is closed with a finite codimension, 

where M u (u,g) denotes the Frechet derivative of M with respect to u. Then, there 
exists sfER and a fi E X 2 that are not both egual to zero such that 

k{Ju(u,g),v) - (n,M u (u,g)v) = 0 VrEli 


and 


min£(«, 2 , pi, k) = C{u, g, pi, k) . 


where C(u, g, pi, k) = kj(u,g) — ( pi,M(u,g )) is the Lagrangian for the constrained 
minimization problem (2.5) and where J u {u,g) denotes the Frechet derivative of J 
with respect to u. Moreover, if 

(2.9) the algebraic sum M u (u, g)X\ + M{u,&) contains 0 E X 2 as an interior 

point, 

then we may choose k = l, i.e., there exists a pi E X 2 such that 


(Ju(u,g),v) - (pi,M u (u,g)v) = 0 VuEli 


and 


min£(«, 2 , pi, 1) = C{u, g, pi, 1) . 


Proof: See [21]. □ 

Next, we apply Theorem 2.2 to the optimization problem (2.4). In doing so, we 
will need the following result. 

LEMMA 2.3. Let the spaces X, Y , and Z and operators T and N be defined as m §2.1. 
For v E X, assume that N'(v) maps X into Z . Then, TN'(v) is a compact operator 


5 



from X to X, and therefore <r(— TA'(i>)) , the spectrum of the operator (— TA'(i>)), is 
at most countable with zero being the only possible limit point. 

Proof: Since Z Y , we see that N'(v) is a compact linear operator from X to 

Y . Also, T is a bounded linear operator from Y to A so that TN'(v) is a compact 
operator from X to X . Hence, <r(— TA'(i>)) is at most countable and consists only of 
0 and the eigenvalues of (-TJV , (o)). □ 

Note that in the following result, the existence of at least one pair (u, g) satisfying 
(2.4) is guaranteed by Theorem 2.1. 

THEOREM 2.4. Let A E a be given. Assume that assumptions (H1)-(H9) hold. Let 
(u,g) E X x 0 be an optimal solution satisfying (2.4). Then, there exists a k E 1R 
and a fi E X* that are not both egual to zero such that 

(2.10) k (J u (u,g),w) - (/j,M u (u,g) -w) = 0 M w E X 
and 

( 2 . 11 ) min £(u, z , /a, k) = £(u, g, /a, k) . 

z£® 

Furthermore, */(l/A) a{—TN'(u}), we may choose k = 1, i.e., for almost all X, 

there exists a fi E X* such that 

(2.12) (Ju(u,g), w) - (fj.,M u (X,u,g) ■ w) = 0 M w E X 
and 

(2.13) min£(u, z, pi, 1) = £{u, g, pi, 1) . 


Proof: Let A £ A be given. To show the existence of k and fi such that (2.10) and 
(2.11) are valid, we only need to verify that the hypotheses (2.6)-(2.8) of Theorem 2.2 
hold with X\ = A 2 = X, since in this case (2.5) reduces to (2.4). Obviously, (2.6) is 
merely a restatement of (H7). Since O is convex and since the mappings T and K are 
linear, we have that if z 7 = 7 Z 1 + (1 — j)z 2 , then 

M (v, z 7 ) = v + A TN(v) + A TKzj 

= j(v + XTN(v)) + (1 - y)(v + XTN(v)) + jX(TKzi + (1 - j)TKz 2 )) 

= 7 M(v, zi) + (1 - j)M(v, z 2 ) . 

Moreover, (H8) implies that 

J{ v, z 7 ) = Af(ti) + A £(z 7 ) = Af(c) + \£(jzi + (1 - j)z 2 ) 

< Af(c) + A (y£(zi) + (1 - 7 )£(z 2 )) = 7 J(v, zi) + (1 - 7 )J(v, z 2 ) . 

Thus, (2.7) holds. The operator M u (u,g) from A to A is defined by 

M u (u, g) ■ w = w + XTN' (u) ■ w M w E X 


or simply, 


M u (u,g) = I + XTN'(u). 
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From (H9) and Lemma 2.3, we have that TN'(u) is a compact operator from X to X . 
As a result, M u (u,g) = I + XTN'(u) is a Fredholm operator so that it has a closed 
range with a finite codimension, i.e. , (2.8) holds. Thus, by Theorem 2.2, there exists 
a k G IR and a fi G X* which are not both equal to zero such that (2.10) and (2.11) 
hold. 

To show the existence of a fi such that (2.12) and (2.13) are valid, we only need 
to verify that the additional hypothesis (2.9) of Theorem 2.2 holds. In fact, if, in 
addition (1/A) ^ <r(— TN'(u)') , then it follows that A = Range(/ + \TN'(u)') = 
Range (M u (u, < 7 )) so that Range(M u («, < 7 )) contains 0 £ A as an interior point, i.e., 
(2.9) holds. Hence, by Theorem 2.2 and Lemma 2.3, we conclude that for almost all 
A, there exists a fi £ X* such that (2.12) and (2.13) hold. □ 

So far 0 has only been assumed to be a closed and convex subset of G. No 
smoothness condition on the control variable g has been assumed in the functional 
or in the constraint. Thus, the necessary condition of optimality with respect to 
variations in the control variable is expressed in the cumbersome relation (2.11). We 
now turn to the case where 0 contains a neighborhood of g, where (u, g) is an optimal 
solution. In particular, we assume that 0 = G. In this case, (2.11) can be given a 
more concrete structure. 

THEOREM 2.5. Let A £ A be given. Assume that assumptions (H1)-(H10) hold. Let 
(u, g) G X x G be a solution of the problem (2.4). Then, there exists a k G IR and a 
fi G X* that are not both equal to zero such that 

(2.14) k {Ju(u, g),w) - {n, (/ + \TN'(u))w) = 0 Vw£A 
and 

(2.15) k{8'{g),z) - {n,TKz) = 0 VzGG. 

Furthermore, */(l/A) ^ <r(— TN'(u)), we may choose k = 1, i.e., there exists a qi G X* 
such that 

(2.16) (Ju(u, g),w) - (j!, (L + \TN'(u))w) = 0 Vw£A 
and 

(2.17) (S'(g),z) - (n,TKz) = 0 VzGG. 
hold. 

Proof: Since the hypotheses imply that J(v,z) is Frechet differentiable with respect 
to z, (2. 14)-(2. 17) follow easily from Theorem 2.4. □ 

Remark. If k = 0, then there exists a fi 7 ! 0 such that 

— (a*, M u (u, g)w) = 0 Vi()£l 

so that the optimality system necessarily has infinitely many solutions. In fact, for any 
G G IR , {C n) is a solution whenever fi is a solution. This creates both theoretical and 
numerical difficulties. Thus, it is of great interest to try to eliminate this situation. 
Fortunately, Lemma 2.3 and Theorem 2.4 tell us that we may set £=1^0 for almost 
all values of (1/A), i.e., except for the at most countable set of values in <r(— TA'( m)). 
□ 
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If the control g enters the constraint in a favorable manner, then we may take 
k = 1 even when (1/A) £ <r(— TX'( m )) . Specifically, we invoke one of the assumptions 
(Hll) and (Hll)'. We then have the following result. 

THEOREM 2.6. Assume that the hypotheses of Theorem 2.5 hold. Assume that if 
(1/A) £ <r(— TX'( m)), then either (Hll) or (Hll)' holds. Then, for all A £ A, there 
exists a fi £ X* such that (2.16) and (2.17) hold. 

Proof: Because of Theorem 2.5, we only need to examine the case (1/A) £ <r(— TX'( m)) 
and show that the algebraic sum M u (u, g)X + M(u, G) = X. If (Hll)' holds, the result 
is a direct application of Theorem 2.2. 

If (Hll) holds, let (1/A) be a nonzero eigenvalue of (— TN’(u )). Then, A is 
also an eigenvalue of (— N' (u)*T*^j with a finite dimensional eigenspace having the 
corresponding eigenfunctions {v* C X* as a basis. We claim that {K*T*v * C 
G* is a linearly independent set. To see this, we assume J/fci a i K*T*v* = 0 with 
ay £ 1R ; this expression can be rewritten as K*T* (J27Li a i v *) = 0- Because each v* 
is an eigenvector, we have (/ + A N' (m)*T*) J/fci a i v* = 0. Thus, the assumption 
(Hll) implies that J/fci a i v i = 0- Since {v*}’ j r ! =1 is an eigenbasis, and is therefore 
a linearly independent set, we have each ay = 0. This shows that {K*T*v*}’ j f_ 1 is 
linearly independent set in G* . Hence, we may choose an orthonormal dual basis 
{ z i}T=i C G such that {z iy I<*T*v*) = Sij . 

Now, let w £ X be given. We choose z = yj/fci { w > v i) z i- Then ( w,Vj ) — 
X(TI<z,v*) = (w, Vj) — A(z, K*T*Vj) = (w, v*) ~Ya=i{ w ’ v i) 6 P = 0 for j = l,...,m. 
Thus, by Fredholm alternatives, there exists a unique v £ X that satisfies (/ + 
XTN' (u))v = w — XTKz, or, (/ + XTNfu^v + XTKz = w; thus, we have shown 
that M u (u, g)X + M(u, G) = X. Hence, by Theorem 2.2, there exists a a £ X* such 
that (2.16) and (2.17) hold. □ 

2.5. The optimality system 

Under the assumptions of Theorem 2.6, an optimal state u £ X , an optimal con- 
trol g £ G, and the corresponding Lagrange multiplier fi £ X* satisfy the optimality 
system of equations formed by (2.2), (2.16), and (2.17). From (2.1) we have that 
J u = XT' and J g = XT', where T' denotes the obvious Frechet derivative. Then, 
(2 . 16)-(2 -17) may be rewritten in the form 

(2.18) n + X [N'(u)]*T*n - X T\u) = 0 in X* 
and 

(2.19) £\g) - K*T*pi = 0 in G* . 

For purposes of numerical approximations, it turns out to be convenient to make the 
change of variable / = T* pi. Then, the optimality system (2.2), (2.18), and (2.19) for 
u £ X, g £ G, and / £ Y* takes the form 

(2.20) u + XTN(u) + XT Kg = 0 in X , 

( + AT t [A f, (ti)n-AT*f , («) = 0 inY*, 


£\g) - K*f = 0 in G* . 


( 2 . 21 ) 

and 

( 2 . 22 ) 
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It will also be convenient to invoke an additional simplifying assumption concern- 
ing the dependence of the objective functional on the control. Specifically, we assume 
that (H12) holds. Then, (2.20)-(2.22) can be rewritten as 


(2.23) 

u + A TN(u) + A TKg = 0 in 

(2.24) 

( + AT* [N’(u)]*(, - A T*T'(u) = 0 

and 


(2.25) 

g - EK*£ = 0 in G . 


Remark. Note that the optimality systems, e.g., (2.23)-(2.25), are linear in the adjoint 
variable Also, note that the control g may be eliminated from the optimality system 
(2.23)-(2.25). Indeed, the substitution of (2.25) into (2.23) yields 

(2.26) u + A TN{u) + A TKEK*£ = 0 in A . 

Thus, (2.24) and (2.26) determine the optimal state u and adjoint state subse- 
quently, (2.25) may be used to determine the optimal control g from This observa- 
tion serves to emphasize the important, direct role that the adjoint state plays in the 
determination of the optimal control. □ 

Remark. Given a ^ £ Y* , it is not always possible to evaluate g exactly from (2.25). 
For example, the application of the operator E may involve the solution of a partial 
differential equation. Thus, although it is often convenient to devise algorithms for 
the approximation of optimal control and states based on the simplified optimality 
system (2.24) and (2.26), in some other cases it is best to deal with the full form 
(2.23)-(2.25). Thus, when we consider approximations of optimal controls and states, 
we will deal with the latter. □ 

Remark. In many applications we have that X* = Y . Since these spaces are assumed 
to be reflexive, we also have that Y* = X. In this case, we have that both u and £ 
belong to A. □ 


3. Finite dimensional approximations of the abstract problem 

In this section we define and analyze algorithms for the finite dimensional approx- 
imation of solutions of the optimality system (2.23)-(2.25); an outline of the definitions 
and results of this section is as follows. 

• In §3.1 , we define the finite dimensional approximate problems 
that we consider. 

• In §3.2, a list of assumptions about the approximate problems 
is given. 

• In §3.3, we quote a result of [6] that we will use to analyze 
approximations obtained as solutions of the approximate prob- 
lems defined in §3. 1-3.2. 

• In §3.4, we provide error estimates for the approximation of 
solutions of the optimality system (2.23)-(2.25) 
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3.1. Formulation of finite dimensional approximate problems 

A finite dimensional discretization of the optimality system (2.23)-(2.25) is de- 
fined as follows. First, one chooses families of finite dimensional subspaces X h £ X , 
(Y* ) h £ Y* , and G h £ G. These families are parameterized by a parameter h that 
tends to zero. (For example, this parameter can be chosen to be some measure of the 
grid size in a subdivision of Q into finite elements.) Next, we define approximate oper- 
ators T h : Y ->■ X h , E h : G* G h , and ( T* ) h : X* ->■ {Y* ) h . Of course, one views 
T h , E h , and (T* ) h as approximations to the operators T, E, and T* , respectively. 
Note that (T*) h is not necessarily the same as (T h )* . The former is a discretization 
of an adjoint operator while the later is the adjoint of a discrete operator. 

Once the approximating subspaces and operators have been chosen, an approxi- 
mate problem is defined as follows. We seek u h £ X h , g h £ G h , and £ h £ (Y* ) h such 
that 

(3.1) u h + \T h N(u h ) + \T h I<g h =0 in X h , 

(3.2) f h + X(T*) h [N'(u h )]*f h - \{T*) h X'{u h ) = 0 in {Y*) h , 
and 

(3.3) g h -E h K*C h = 0 in G h . 

3.2. Hypotheses concerning the abstract problem and the approximate 
problem 

We make the following hypotheses concerning the approximate operators T h , 
( T *) h , and E h : 

(H13) lim \\(T — T h )y\\x = 0 VyeY, 

h — ► 0 

(H14) lim \\(T* — (T*) h )v\\y = 0 V v £ X* , 

h — ^0 

and 

(H15) lim 11(1? — E h )s\\a = 0 V s £ G* . 

ft— ► 0 

We also need the following additional hypotheses on the operators appearing in 
the definition of the abstract problem (2.4): 

(H16) N £C 3 (X;Y) and X £C 3 (X;] R); 

(HIT) N" , N'" , T" , and T'" art locally bounded, i.e., they map bounded sets to 

bounded sets; 

(H18) for v £ X, in addition to (H9), i.e., N'(v) £ C(X;Z) where Z Y, 
we have that [tV'(i;)]* £ C(Y*;Z) where Z X* , that for ry £ Y* , 

[At"(i>)]* • r} £ C(Y*\ Z), and that for w £ X , X"{y) ■ w £ C(X; Z); and 
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(H19) K maps G into Z . 

Here, (•)" and (•)'" denote second and third Frechet derivatives, respectively. 

3.3. Quotation of results concerning the approximation of a class of non- 
linear problems 

The error estimate to be derived in Section 3.4 makes use of results of [6] and [10] 
(see also [13]) concerning the approximation of a class of nonlinear problems. These 
results imply that, under certain hypotheses, the error of approximation of solutions 
of certain nonlinear problems is basically the same as the error of approximation of 
solutions of related linear problems. Here, for the sake of completeness, we will state 
the relevant results, specialized to our needs. 

The nonlinear problems considered in [6], [10], and [13] are of the following type. 
For given A £ A, we seek ip £ X such that 

(3.4) n(\,ip) = ip + Tg(\,ip) = o, 

where T £ £(JA;T), Q is a C 2 mapping from Axf into y, X and y are Banach 
spaces, and A is a compact interval of 1R . We say that {(A, -0(A)) : A £ A} is a branch 
of solutions of (3.4) if A — s- 0(A) is a continuous function from A into X such that 
7i(X, 0(A)) = 0. The branch is called a regular branch if we also have that 70 (A, 0A)) 
is an isomorphism from X into X for all A £ A. Here, ?0(-, •) denotes the Frechet 
derivative of %{■, •) with respect to the second argument. We assume that there exists 
another Banach space Z, contained in y, with continuous imbedding, such that 

(3.5) g^ip) e C(X;Z) V A £ A and if £ X, 

where £0(-, •) denotes the Frechet derivative of (/(•, •) with respect to the second ar- 
gument. 

Approximations are defined by introducing a subspace X h C X and an approxi- 
mating operator T h £ T(T; X h ). Then, given A £ A, we seek ip h £ X h such that 

(3.6) n h (\,ip h ) = ip h + T h g(\,ip h ) = o. 

Concerning the operator T h , we assume the approximation properties 

(3.7) lim \\(T h — T)lo\\x = 0 Vw£T 

h — ^0 

and 

(3.8) hrn\\(T h -T)\\ c(Z;X) = 0. 

Note that whenever the imbedding Z C T is compact, (3.8) follows from (3.7) and, 
moreover, (3.5) implies that the operator TC0( A, ip) £ C(X] X) is compact. 

We can now state the result of [6] or [10] that will be used in the sequel. In the 
statement of the theorem, D 2 g represents any and all second Frechet derivatives of g . 

THEOREM 3.1. Let X and y be Banach spaces and A a compact subset of 1R . Assume 
that g is a C 2 mapping from A x X into y and that D 2 g is bounded on all bounded 
sets of A x X . Assume that (3.5), (3.7), and (3.8) hold and that {(A, 00); A £ A} is a 
branch of regular solutions of (3.4). Then, there exists a neighborhood O of the origin 
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in X and, for h < ho small enough, a unique C 2 function A — >■ t/> ft (A) G such that 
{(A, t/> ft (A)); A G A} is a branch of regular solutions of (3.6) and (A) — t/>(A) G O for 
all A G A. Moreover, there exists a constant C > 0, independent of h and A, such that 

(3.9) || t /’ ft (A) — "0(A) 1 1 a* < C'\\(T h — T)Q(\, t/’(' S 0)IU VAgA. □ 

3.4. Error estimates for the approximation of solutions of the optimality 
system 

We now apply the result of Theorem 3.1 to study the approximation of solutions 
of the optimality system. Set X = XxGxY*,y = Yx X* , Z = Z x Z , and 
X h = X h x G h x (Y*) h . (Recall that Z was introduced in (H18).) By the hypotheses 
on Z and Z, we have that Z is compactly imbedded into y. Let T G £(T;T) be 
defined in the following manner: T(r, f) = ( u,g , £) for (r , f) G T and (u,g,£) G X if 
and only if 

(3.10) u + Tr = 0, 

(3.11) i+T*f = 0, 
and 

(3.12) g-EK*i=0. 

Similarly, the operator T h G C(y\X h ) is defined as follows: T h (r,r) = (u h ,g h ,£ i h ) 
for (r , f) G y and (u h , g h , ^ h ) G X h if and only if 

(3.13) u h + T h r = 0, 

(3.14) i h + (T*) h f = 0, 
and 

(3.15) g h -E h K*i h = 0. 

The nonlinear mapping Q : A x X — >■ y is dehned as follows: Q (A, (w, g, £)) = (r , f) 
for A G A, (m, j,() G and (r , f) G T if and only if 

(3.16) r = XN(u) + XKg 
and 

(3.17) f = A[W(h)]*f- AJ^h). 

It is evident that the optimality system (2.23)-(2.25) and its finite dimensional coun- 
terpart (3.1)-(3.3) can be written as 

(u,g,0 + T G( x A u ,g,0) = 0 


12 



and 


(u h , 9 h ,e)+T h g(\,(u h ,g h ,e)) = o, 

respectively, i.e. , with ip = (u,g,f) and if) h = (u h ,g h in the form of (3.4) and 
(3.6), respectively. 

Now we examine the approximation properties of T h . 

LEMMA 3.2. Let the operators T and T h be defined by (3.10)-(3.12) and (3.13)-(3.15), 
respectively. Assume that the hypotheses (H13)-(H15) hold. Then, 

(3.18) lim ||(T — T h )(r, t)\\x = 0 V (r, r) E y . 

ft— ► 0 


Proof: Let (u,g,£) = T(r,r), i.e., ( u,g ,£ ) satisfies (3 . 10)-(3 . 12) . Let (u h ,g h ,£ h ) = 
T h (r,r) , i.e., ( u h ,g h ,£ h ) satisfies (3 . 13)-(3 . 15) . Subtracting the corresponding equa- 
tions yields that 

\\u-u h \\ x = \\(T-T h )r\\ x , 

\\i-i h \\Y* = \\(T*-(T*) h )T\\ Y *, 

and 

\\9 ~ 9 h \\a = || (E - E h )K*£ h + EL<*(£ - | h )|| G 

< || (E - E h )K*£ h \\ G + \\EK*\\ £(y , ]G) ||(f - ?)\\ G . 

Thus, for some constant C > 0, 

\\(T-T h )(r,r)\\ x 

< C {||(T - T h )r\\ x + || (T* - (T*) h )r\\ Y , + \\(E - E h )L<*£ h \\ G } . 

Then, the result of the proposition follows from (H13)-(H15). □ 

Next, we examine the derivative of the mapping Q . 

LEMMA 3.3. Let the mapping Q : Ax X — >■ y be defined by (3 . 16) -(3.17). Assume that 
the hypotheses (H9), (H16), and (H18)-(H19) hold. Then, for every A £ A and every 
(u,g,f ) E X, the operator £(«,<,,£) (A, («, ff, 0) e C(X]Z). 

Proof: A simple calculation shows that G( u ,g,()( A, (u,g, 0) ^ £( x ',y) is gi yen % 

Q(u , 9 ,0 (A, (u, g, 0) • («, g, £) = A ( . „]* . £ | [N'(u)]* A? - E"{u) ■ u ) ' 

Then, the result follows from (H9) and (H18)-(H19). □ 

A solution (u(A), g(X), £(A)) of the optimality system (2.23)-(2.25) is called regular 
if the system (for the unknowns (u,g,£)) 

(3.19) u + ATAf'(«)fi + XTKg = x, 

(3.20) £ + A T* [#"(«)]*« - ( + AT* [N’(u)]*£ - A T*T"(u)u = y 
and 

(3.21) g — EL{*£ = z . 
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is uniquely solvable for any ( x,z,y ) £ X = X x G x Y* . (Note that the linear 
operator appearing on the left hand side of (3.19)-(3.21) is obtained by linearizing the 
optimality system (2.23)-(2.25) about (u,g,f).) 

In the following theorem, we will assume that the solution («( A), 17 (A), £(A)) of the 
optimality system (2.23)-(2.25) that we are trying to approximate is a regular solution. 
The assumptions we have made, in particular (H9), (H18)-(H19), are sufficient to 
guarantee that for almost all values of A, this is indeed the case. 

LEMMA 3.4. Assume ihe hypotheses of Lemma 3.3. Then, for almost all X, solutions 
(«(A), 77 (A), £(A)) of the optimality system (2.23)-(2.25) are regular. 

Proof: The system (3.19)-(3.21) is equivalent to 

(3.22) (/ + XTS(u,g,f)^j(u,g,£) = ( x,z,y ), 
where the linear operator S(u,g,f) : X —>■ y is defined by 

•S(u,g,f) ■ (u, < 7 , £) = j0(u, g ,o( X ’ ( u ,9,0) ' 

_ f N'(u) ■ u + Kg 

~ ^ [#"(«) •«]*•£ + [N'(u)]* ■ £ - T"(u) ■ u 

Now, T £ T(T; T), so that, by Lemma 3.3, (/ + A T S(u, g , ^)) is a compact perturba- 
tion of the identity operator from X to X . Thus, for almost all A, (3.22), or equivalently 
(3 . 19)-(3 .21), is uniquely solvable, i.e. , for almost all A, the solution («(A), 17 (A), £(A)) 
of the optimality system (2.23)-(2.25) is regular. □ 

Using Theorem 3.1, we can now provide an error estimate for approximations of 
solutions of the abstract problem. 

THEOREM 3.5. Let («( A), 17 (A), £(A)) £ X , for A £ A, be a branch of regular solutions 
of the optimality system (2.23)-(2.25). Assume that the hypotheses (H13)-(H19) hold. 
Then, there exists a 8 > 0 and an ho > 0 such that for h < ho, the discrete optimality 
system (3 . 1) -(3 .3) has a unigue solution (u h (X), g h (X), f h (X)) satisfying 

||(«(A),< 7 (AU(A)) - (u\x),g\x),ew)\\x <s. 

Moreover, 

(3.23) lim || («(A), 77 (A), £(A)) - (« ft (A), g h (X), £ h (A)) ||* = 0 

uniformly in A £ A and there exists a constant C', independent of h and X, such that 

lim || («(A), 77 (A), e(A)) - (« ft (A),/(A),^(A))|U 

h — ► 0 

(3.24) < UA{||(T ft - T)(N(u( A)) + A'^A)) || x + \\(E h - E)L<*f(X)\\ G 

+ ll((^*) ft - T*) {[N\u(X))n ~ ^'(«(A))) ||y* } • 



Proof: Assumptions (H16) and (H17) ensure that Q £ C 2 (X,y) and D 2 Q maps 
bounded sets of A x X into bounded sets of y . By Lemma 3.3, assumptions (H18) 
and (H19) imply that (3.5) holds. By Lemma 3.2, assumptions (H13)-(H15) imply 
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that (3.7) holds. Then, since Z is compactly imbedded into y , (3.7) implies that (3.8) 
holds. Thus, all the hypotheses of Theorem 3.1 are verified. Then, a direct application 
of Theorem 3.1 yields (3.23) and (3.24) follows from (3.9). □ 

It is easily seen that (3.23) and (3.24) are equivalent to: 

lim{|KA) - « ft (A)|| x + || fl (A) - g\X)\\ G + ||£(A) - £ h (A)||y. } = 0 

uniformly in A 6 A and that there exists a constant C , independent of h and A, such 
that 


|KA) - u h ( A)|| x + || fl (A) - g h (X)\\ G + ||£(A) - £ h (A)||y. 

< Cx{\\(T h - T)(N(u( A)) + Kg(X))\\x + \\(E h - A)A^(A)|| G 

+ ii (cn h - t*) ([N'(u(x))rax) - ^w*))) u y . } . 

If, in (3.9), the operator T is invertible, we have, using (3.4), that 
||^(A) - iP(X)\\ x < C\\(T h l 1 - I)iP(X)\\ x V A G A . 

Thus, if the operator T from Y to X is invertible, we have that (3.24) is equivalent to 
|KA) - « ft (A)|| x + || fl (A) - g\X)\\ G + ||£(A) - € h (A)||y. 

(3.25) < C'jlKT'T- 1 - I)u{ A)|| x + | K^if 1 - /) fl (A)|| G 

+ ||((T*) ft (T*)- 1 -7)e(A)|| y .}. 


4. Applications 


We now apply the framework and analyses developed in §2 and §3 to some concrete 
problems, all of which feature constraints on admissible states and controls that take 
the form of a system of nonlinear partial differential equations. In each application, 
we use a different control mechanism so that the discussion provided in this section 
illustrates the treatment of a variety of such mechanisms. However, one could use any 
of the control mechanisms discussed in any of the applications in any other application, 
or in fact, use any combination of such mechanisms. 

Before examining any specific application, we establish some notation. Further 
notation will be established as needed when the individual applications are considered. 

Throughout, C will denote a positive constant whose meaning and value changes 
with context. Also, H 4 S (T>) for s £ R denotes the standard real Sobolev space of order 
s with respect to the set V, where V could either be a bounded domain £2 E IR d , 
d = 2,3, or part of the boundary T of such a domain. Of particular interest are the 
spaces H°(V) = L 2 (V), 


H 1 (V) = \^eL 2 (V) | ^eL 2 (V) for j = 1, d 


and 


*»«>>- | **m-i <>■ 
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Also of interest is the subspace 

H^V) = {<j> e H\V) 


= 0 on 


dV } 


where dV denotes the boundary of V . 

Dual spaces will be denoted by (•)*. Duality pairings between spaces and their 
duals will be denoted by (•, •}. Norms of functions belonging to H s (£l) and H S (T) are 
denoted by || • || s and || • || S) r, respectively. Of particular interest are the F 2 (£l)-norm 
|| • 1 1 o , the f? 1 (Sl)- norm 


and the f? 2 (Sl)-norm 


IHI? = E 

j = i 


Ml= E 

j,k = l 


d „ ^ ll 2 


dx; 


d 2 


\mi 


dxjdxj. 


IHI?- 


Corresponding Sobolev spaces of real, vector-valued functions having r compo- 
nents will be denoted by H S (T>), e.g., H 1 (X>) = [H 1 (V)] r . Of particular interest will 
be the spaces L 2 (£>) = H°(£>) = [L 2 (V)] r , 

f c)v ■ 

H 1 (T>) = < Vj £ L 2 (V) — — £ L 2 (V) for j = 1, . . . , r and k = 1, . . . , d 

l ox k 


and 


n 2 (v) = \ Vj £ l 2 (v) 


92 Vj 

dx k ’ dx k dxi 


£ L 2 (n) 


for j = 1, ... ,r and k,£ = 1 , ,d 


where Vj , j = 1, . . . , r, denote the components of v. Also of interest is the subspace 

Hj(£>) = | v E H 1 (X>) Vj = 0 on dV , j = 1, . . ,,r } . 

Norms for spaces of vector-valued functions will be denoted by the same notation as 
that used for their scalar counterparts. For example, 


= Ell^'H* and H v ll?,r = El 


O lU.r ■ 


i = i 


i = i 


We denote the L 2 (Q.) and L 2 (S1) inner products by (•, •), i.e. , for p,q £ L 2 (Q.) and 


u, v £ L 2 (0) 


(p, q) = / pqdi2 and (u,v)= / u • v di2 . 

J Cl J Cl 


Similarly, we denote by (•, -)r the L 2 (V) and L 2 (T) inner products, i.e., for p,q £ L 2 (V) 
and u, v £ L 2 (T) 

(p, g) r = J pqdT and (u,v)r = J u • v dr . 
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Since in all cases L 2 -spaces will be used as pivot spaces, the above inner product 
notation can also be used to denote duality pairings between functions defined on 
H s - spaces and their dual spaces. 

For details concerning the notation employed, one may consult, e.g., [1], 

4.1. Distributed controls for the von Karman plate equations 


For this application we will use distributed controls, i.e. , control is effected through 
a source term in the governing partial differential equations. Let SI be a bounded, con- 
vex polygonal domain in 1R 2 and let T denote the boundary of SI. The von Karman 
equations for a clamped plate are given by (see, e.g., [9] or [18]) 


A 2 1pl + —[-02, -02] = 0 


in SI , 


and 


where 


A 2 0 2 — [01, 02] = Agr in SI , 


/ i d ^ 2 n r 

Wi = -x— = W 2 = -x— = 0 on r , 

on on 


d 2 i l> d 2 <f> <9 2 i l> d 2 <f> d 2 i l> <9 


r i ^ y ^ r . ^ y 

^ “ dx 2 dxl + dxl 


2 dx l 


- 2 


< 9 * 1*2 < 9 * 1*2 


Here, i/>i denotes the Airy stress function, 02 the deflection of the plate in the direction 
normal to the plate, Xg is an external load normal to the plate which depends on the 
loading parameter A, and <9(-)/<9 n the normal derivative in the direction of the outer 
normal to T. 

By introducing appropriate rescalings, i.e., by replacing 0i by A0i , 02 by A02, 
and g by A g, we can rewrite the von Karman equations as follows: 


(4.1) 

A 2 0 + -^-[02, 02] = 0 in SI, 

(4.2) 

A 2 i />2 - A [0i , 0] = \g in SI 

and 


(4.3) 

, 001 , 002 „ 

01 = = 1 p 2 = - 7 T- = 0 on 

on on 


We introduce the spaces 

H 2 (n) = (0 e h 2 ( n) 


dih 

0 = 0, — = 0 on r 
on 


H 2 (S1) = [F 2 (S1)] 2 , H- 2 (Cl) = (Hq({ 2))* , and H- 2 (S1) = (H 2 (fi))* 

and the bilinear form 


j(0, 0) = f Ai/> A<f> dSl V i>,<t> E tf 2 (Sl) 

Jn 
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in order to define the following weak formulation of the von Karman equations (4.1)- 
(4.3): find ip = (t/q, 1 P 2 ) £ H 2 (ff) such that 

(4-4) a ( t /’i; <Pi) + (["02 , V’ 2 ], <Pi) = 0 V <£1 £ -ffo(fi) 

and 

(4-5) a(-02, <^ 2 ) — A([r/>i, -02], <^ 2 ) = A(gr, <^ 2 ) V </> 2 £ H 2 (£l) . 

Using the identity 

(4.6) (XX\X) = (bP, CU) ^ P’XX £ Hq(Q) , 

one can show that for each g £ H~ 2 (£l), (4.4)-(4.5) possesses at least one solution 
ip = (t/q, ip 2 ) £ H 2 (ff) and that all solutions of (4.4)-(4.5) satisfy the a priori estimate 

(4.7) IXih + WXh <C\\g\\- 2 ; 

see, e.g., [18], for details. In the sequel a solution to (4.1)-(4.3) will be understood in 
the sense of (4.4)-(4.5). 

Given a desired state ip 0 = (ipio, V’ 20 ) £ L 2 (f2), we define for any ip = (t/q, t/q) £ 
Hg(n) and g £ L 2 (Q) the functional 

JA,g) = JAiA2,g) 

^ 8 ' ) / ((Vh - V’lo) 2 + (i >2 - V’ 20 ) 2 ) dtt + ^ 

We then consider the following optimal control problem associated with the von 
Karman plate equations: 

(4.9) min{ J(ip,g) \ ip £ H 2 (ff) , g £ 0 } subject to (4.4)-(4.5) , 
where 0 is a subset of L 2 (X). 

We define the spaces X = H 2 (£l), Y = H _2 (£l), G = L 2 (QP), and Z = L 1 (S1). 
By compact imbedding results, Z Y . For the time being, we assume that the 

admissible set 0 for the control g is a closed, convex subset of G = L 2 (X). 

Let the continuous linear operator T £ C(Y]X) be defined as follows: for f = 
(/ 1 , f 2 ) £ Y = H _2 (S1), Ti = ip £ X = H 2 (f!) is the unique solution of 

a(XiXi) = (fiXi) V pi £ Hl{QP) 


and 

a(X, h) = (/2, fa) V«/> 2 £f7 2 (fi). 

It can be easily verified that T is self-adjoint. 

We define the (differentiable) nonlinear mapping N : X —>■ Y by 


N(ip) = 


\bP2A2] 

-[■P1A2] 


Vipex 
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or equivalently 


(N(tl>),<j>) = -([V’2,'02],^l) - ([lpl,1p2],fa) V ^ = (<0,<?C) G X 


and define K : g £ T 2 (ft) — >■ Y by 


Kg = - 


or equivalently, 


(Kg,4) = 2 } v^ = (0i, 0 2 ) ex. 

Clearly, the constraint equations (4.4)-(4.5) can be expressed as 

ip + X TN(ip) + X TKg = 0, 

i.e. , in the form (2.2). With the obvious definitions for tF(-) and £(•), i.e. , 
Hi>) = \J ((^ - V’io) 2 + (1P2 - ^2o) 2 ) dn vj>£X 

and 

£(a) = \ [ a 2 dn VgtG, 

the functional (4.8) can be expressed as 

J(iP,g) = XX(iP) + X£(g), 


i.e., in the form (2.1). Thus, the minimization problem (4.9) is in the form of the 
minimization problem (2.3). 

We are now in a position to verify, for the minimization problem (4.9), all the 
hypotheses of §2 and §3. 


4-1.1. Verification of the hypotheses for the existence of optimal solutions. 

We first verify that the hypotheses (H1)-(H6) hold in the current setting. 

(HI) is obviously satisfied with a lower bound 0. 

(H2) holds with a = 1 and j3 = 2. 

(H3) is verified with the choice g*- 0 - 1 ) £1x0, where g 1 - 0 ) is an arbitrarily 

chosen element in 0 and ip^ = ('f/’i°\ '02°'') i s a solution of 

4^ 2 0i O) + ■^-[■02° ) . ■*/ , 2° ) ] = 0 infi ) 


A 2 ip 


( 0 ) 

2 




4 0> 


XgG) in ft, 


and 


( 0 ) _ dip 


( 0 ) 


tr = 


dri 


( o) _ dipj 0) _ 


00 = 


dri 


= 0 


on r . 


In order to verify (H4), we assume C 0 is a sequence satisfying g( n ) g 

in T 2 (ft); then, we have gG) g in H~ 2 (Vt ) so that lim n ^ 00 (^r( n ), z) = ( g,z ) for 
all z £ H 2 (Vt), i.e., KgG) Kg in Y. Assume that the sequence {ip n ’ > } C H 2 (ft) 


19 



satisfies 0 n - ) if) in H 2 (fl); then, (d 2 ij) ( '' n ' j / dxidxj) (d 2 ij) / dxidxj) in L 2 (£!) and, 

by using a compact imbedding result, 0 n - ) — >■ in L 2 (S1). Now, using the identity 

(4.6), 

(N(4> (n) ),<t>) = \ ([4 n) ,4 n) ]^i) - 

= \ - ([V’i n) ,^2],t/’2 n) ) 

^ ([0, 4>i], 0) - ([0, fa], 0) 

= ^ ([ 0 , 0 ], 0 ) - ([ 0 , 0 ], 0 ) = . 

Hence, (H4) is verified. 

The verification of (H5) follows directly from the observation that the mappings 
4 ^ F(4) = M/MW4 - 0||o and g ^ £(g) = (l/ 2 )|| 5 f||g are convex. 

The verification of (H6) is a trivial consequence of the a priori estimate (4.7). 

It is now just a matter of citing Theorem 2.1 to prove the existence of an optimal 
solution that minimizes (4.8) subject to (4.4)-(4.5). 

THEOREM 4.1. There exists a (<f>, g) £ H 2 (fl) x 0 such that (4.8) is minimized subject 
to (4.4)-(4.5). □ 


4-1.2. Verification of the hypotheses for the existence of Lagrange multipliers. 

We now assume (ij>, g) is an optimal solution and turn to the verification of hy- 
potheses (H7)-(H9). 

The validity of (H7) is obvious. 

(H8) holds since the mapping g 1 — > £{g) = (1/2)||^|| 2 is convex. 

(H9) can be verified as follows. For any if) £ X , the operator N'(fi)) : X —>■ Y is 
given by 




[■02, 02] 

[ 01 , 02 ] - [ 02 , 0 ] 


V 4 — (0i , <t> 2 ) £ x . 


Thus, using the definition of [•, •], we obtain that N'(tf) ■ $ £ L 1 (S2) = Z . 
The Lagrangian is given by 


£( 0 ,fib»? 0 ) = - <j MM 1,771) + — ([02, 02] , 171) 


MM 2,m) - ^([0i, 02 ], m) - Mg, m)\ 


for all (0 g, r),k) £ X x G x X x 1R = H 2 (fi) x L 2 (tt) x H 2 (fi) x 1R . Note that in this 
form of the Lagrangian, the Lagrange multiplier r) £ X = Y* so that we have already 
introduced the change of variables indicated between (2.17 )-(2 .18) and (2 . 19)-(2 .21). 

Having verified the hypotheses (H7)-(H9), we may apply Theorem 2.4 to conclude 
that there exists a Lagrange multiplier r) £ X = H 2 (fl) and a real number k such that 

(4.10) r) + XT* ([A0VOr 0) = 0 

and 

(4.11) C{4,g,r),k) < C{4,z,r),k) Vz £ 0 
and that for almost all values of A, we may choose k = 1. 
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Recall that T is self-adjoint. Also, note that for any if) £ X = H 2 (£2), 




-[7/2,772] 
[7/2,771] - [7/1,772] 


v 17 = (771,772) ex . 


Thus, (4.10), with k = 1, can be rewritten as 

(4.12) a(C 1 ,77 1 )-A([^ 2 , 77 2 ],C 1 )=A(^ 1 -^ 10 ,Ci) V (i G H 2 0 (n) 

and 

(4 13 . a ((2 , 772) + A ([V>2 , 771] , (2) - A ([ipl , 772] , C2) 

= A(l/ , 2 7 ^ 20 , C2) V ( 2 eH 2 ( n ). 


Using the definition of the Lagrangian functional, (4.11), with k = 1, can be rewritten 
as 

- (z, z) + A(z, 772 ) - - (g, g) - \{g, t? 2 ) > 0 V z E 0 . 

Note that, in the above expression, we have already employed hypothesis (H12) which 
in the current context is trivially satisfied with E the identity operator on G* = G = 
L 2 (£2). For each e £ (0, 1) and each t £ 0, set z = et + (l — e)g £ 0 in the last equation 
to obtain 

e 2 

y (t ~ 9,t - g) + e (t - g, g) + e (t - g, m) > 0 V t £ 0 

so that, after dividing by e > 0 and then letting e —>■ 0 + , we obtain 
(4.14) (t - g,g + t? 2 ) > 0 Vt£0. 


We see that for almost all values of A, necessary conditions for an optimum are 
that (4.4)-(4.5) and (4.12)-(4.14) are satisfied. The system formed by these equations 
will be called an optimality system. 

We now specialize to the case 0 = T 2 (£2). Note that the hypothesis (H10) is 
satisfied. Then, using Theorem 2.5, we see that the inequality (4.14) becomes an 
equality and, by letting z = t — g vary arbitrarily in T 2 (S2), we now have, instead of 

(4.14) , 

(4.15) (z, g + 772 ) = 0 V z £ f 2 ( f2) . 


Thus, according to that theorem, we have that for almost all A, an optimality system 
of equations is now given by (4.4)-(4.5), (4. 12)-(4. 13) , and (4.15). However, we can go 
further and verify that the hypothesis (Hll)' is valid, which in turn will justify the 
existence of a Lagrange multiplier satisfying the optimality system for all A £ A. We 
now assume the domain £2 is a convex polygon with no angles greater than 126°. 

Let A be given such that 1/A is an eigenvalue of —TN'(ij>), where ( if),g ) £ H 2 (£2) x 
L 2 (£2) is an optimal pair that minimizes (4.8) subject to (4.4)-(4.5). We wish to show 
that for each f £ H -2 (£2), there exists a g £ L 2 (Q) and a i/> £ H 2 (£2) such that 


if) + \ TN'(if)) ■ if) + A TKg = f , 


i.e., 

( 4 . 16 ) 0(7/1, <f>i) + A([V> 2 , 7/2], <f>i) = (/i,</i) V cf>i £ f?o(£ 2 ) 
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and 


(4.17) 


a(f> 2, fa) ~ A ([V>i, V> 2 ], <^2) - A([V>i, V> 2 ], ^2) - A(fif, <^2) 

= (/2,^2) v <p 2 e H 2 (n) . 


To show this, we first let ip £ Hq(S 1) be a solution of 


a(i>i,4>i) + A([V>2,'0 2 ],^i) = (fiAi) V (f>! £ H%(Q) 


and 

a(i> 2 , fa) - a([V>i, V> 2 ], </> 2 ) = {f 2 , fa) v^e H 2 (n) . 

The existence of such a ip can be shown in a manner similar to that for showing the 
existence of a solution to the von Karman equation; the key step is that by adding the 
two equations with the test function <f> replaced by ip, we have the a priori estimate 

+ a (i , 2, V> 2 ) = (/l,'0l) + (/2,'0 2 ) • 

Then, we choose g = — [t/q, ip 2 ]- Note that regularity results for the biharmonic equa- 
tion applied to (4.4)-(4.5) yield rp £ H 4 (fl) (see [3]). Hence, using imbedding theorems 
we deduce that g £ L 2 ( SI). It is obvious that g and ip satisfy (4. 16)-(4. 17) , i.e. , we have 
verihed (Hll)'. Hence we conclude that for all A, the optimality system (4.4)-(4.5), 
(4. 12)-(4. 13) , and (4.15) has a solution. Thus, we have Theorem 2.6 which, in the 
present context, is given as follows. 

THEOREM 4.2. Let (ip, g) £ Hq(S2) x L 2 (Ll) denote an optimal solution that minimizes 
(4.8) subject to (4.5)-(4.6). Then, for all X £ A, there exists a nonzero Lagrange 
multiplier r) £ Hq(S2) satisfying the Euler equations (4.12)-(4.13) and (4.15). □ 

4.1.3. Verification of the hypotheses for approximations and error estimates. 

We finally verify the hypotheses (H13)-(H19) that are used in connection with 
approximations and error estimates. 

A finite element discretization of the optimality system (4.4)-(4.5), (4. 12)-(4. 13) , 
and (4.15) is defined in the usual manner. We first choose families of finite dimensional 
subspaces X h C Hq(S 2) and G h C L 2 (Ll) parameterized by a parameter h that tends 
to zero and satisfying the following approximation properties: there exists a constant 
C and an integer r such that 

(4.18) inf \\<P-4> h \\ 2 <Ch m \\4>\\ m+2 , V<f> £ H m+2 (fi), 1 < m < r 
<f> h £X h 

and 

(4.19) inf \\z-z h \\ 0 < Ch m \\z\\ m , V z e H m (Q ) , 1 < m < r . 

z h £G h 

One may consult, e.g., [8] for some finite element spaces satisfying (4.18) and (4.19). 
For example, one may choose X h = V h x V h where V h is the piecewise quintic-C ,1 (S2) 
finite element space constrained to satisfy the given boundary conditions and defined 
with respect to a family of triangulations of S2. In this case, h is a measure of the grid 
size. For simplicity, one may choose G h = V h . 
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Once the approximating spaces have been chosen, we may formulate the approx- 
imate problem for the optimality system (4.4)-(4.5), (4.12)-(4.13), and (4.15): seek 
ij> h £ X h , g h £ G h , and T) h £ X h such that 

(4.20) a(^i^i)+f([^2.^2]^i) =0 

(4.21) a(^,^)-A([^,^],^) = ( fl ,^) V h , 

(4.22) a(Ctni)- \([iP h 2 ,riUi)=Mrt-Tpio,(i) V Ci G V h , 

(4.23) a(C 2 h ,7, 2 h ) + A([^,^],C 2 h ) -A([^,ry 2 ft ],C 2 ft ) = \(^ 2 ~1>20,$) ^(2 e V h , 

and 

(4.24) (z h ,g h +rfo) = 0 V z h £ G h . 

The operator T h £ C(Y;X h ) is defined as follows: for f £ Y, T h f = tf) h £ X h is 
the unique solution of 


a(4> h l A h i) = (hA h i) y^ev h 

and 

«(^ 2 ^ 2 ) = (/ 2 ^ 2 ) y^eV k . 

Since T = T* , we define (T* ) h = T h . 

We dehne the operator E h : L 2 (Q) — >■ G h as the T 2 (£2)-projection on G h , i.e. , for 
each g £ L 2 (Q), 

(E h g^ h ) = (g^ h ) y ij> h £ G h . 

Since G = L 2 (£l) is reflexive, E h is in fact an operator from G* — >■ G h . 

By the well-known results concerning the approximation of biharmonic equations 
(see, e.g., [2] or [8]), we obtain 


||(T — T h )f\\x — >■ 0 

as h —>■ 0, for all f £ Y . This is simply a restatement of (H13). 

(H14) follows trivially from (H13) and the fact that T is self-adjoint and we have 
chosen (T*) h = T h . 

(H15) follows from the best approximation property of T 2 (£2)-projections and 
(4.19). 

(H16) and (H17) follow from the fact that N and T are polynomials. Here we 
also use imbedding theorems and Cauchy inequalities. 

We set Z = Z = L 1 (S2). For each r) £ H 2 (£2) and £ £ H 2 (£2), Sobolev imbedding 
theorems imply that 




~[^2,m] \ a nr 

[^2,m] - 
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and 

■« ■’>=(*<!) € ^' 

These relations verify (H18). 

From the definition of the operator K we see that K maps L 2 (Q.) into 1/(12), i.e. , 
K maps G into Z . Thus (H19) is verified. 

Hence, we are now in a position to apply Theorem 3.5 to derive error estimates 
for the approximate solutions of the optimality system (4.4)-(4.5), (4.12)-(4.13) and 
(4.15). It should be noted that Lemma 3.4 implies that for almost all values of A, the 
solutions of the optimality system are regular. 

THEOREM 4.3. Assume that A is a compact interval of 1R_|_ and that there exists 
a branch {(A, ^>(A), g(A), 17 (A)) : X £ A} of regular solutions of the optimality system 
(4.4)-(4.5), (4.12)-(4.13), and (4.15). Assume that the finite element spaces X h and 
G h satisfy the hypotheses (4.18)-(4.19). Then, there exists a 6 > 0 and an ho > 0 such 
that for h < ho, the discrete optimality system (4.20)-(4.24) has a unique branch of 
solutions {(A, xj} h (X), g h (X), ij ft (A)) : A £ A} satisfying 

m h W - V>(A )|| 2 + H/(A) - <?(A) 1 1 o + ll^(A) - 17 (A) 1 1 2 } < ^ for all X £ A. 

Moreover, 

lim{||^(A) - V>(A )|| 2 + \\g h (X) - (A) 1 1 0 + | |i? ft (A) - »?(A)|| 2 } = 0 , 

h — ► 0 

uniformly in A £ A. 

If, in addition, the solution of the optimality system satisfies (^>(A), g(X), 17 (A)) £ 
H m+2 (f2) x H m (Tf) x H m+2 (f2) for X £ A, then there exists a constant C, independent 
of h, such that 


IIV’(A) - i> h (X)\\ 2 + |b(A) - g h (X)\\ 0 + ||i?(A) - i? ft (A )|| 2 
< Ch m (||V>(A )|| m+2 + ||ff(A)|| m + ||i?(A)|| m+2 ) , 


uniformly m A £ A. 

Proof: All results follow from Theorem 3.5. For the last result, we also use (3.25) and 
the estimates (see, e.g., [ 2 ] or [ 8 ]) 

IIC^T - 1 - 7)V >|| 2 < Ch m \\i>\\m +2 for V e H m+2 (fl) , 

||((T*) h (T*)-i - I)n || 2 = IKT^T - 1 - I)r,\\ 2 < Ch m \\r,\\ m+2 for r, £ H m + 2 (L1) , 

and 

WiEbE- 1 - I)g\\o < Ch m \\g\\ m for g £ H m (Q ) . 

In these estimates, the constant C is independent of h , if), g, 1 ), and A. □ 

Remark. In fact, we obtain from (4.15) that g = —rj 2 so that the term ||fif(A)|| m in the 
right-hand side of the error estimate is redundant. □ 

Remark. By using (4.15) again, along with (4.24) and the error estimate in Theorem 
4.3, we have the following improved error estimate for the approximation of the control 

9 || fl (A) - /(A )|| 2 = IMA) - 172 (A) 1 1 2 < C'h m (M A )|| m+2 + ||*?(A)|| m + 2 ) • 
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Of course, we also use the fact that we have chosen G h = V h C i7 2 (£2). □ 

4.2. Neumann boundary controls for the Ginzburg-Landau superconduc- 
tivity equations 

For this application we will use Neumann boundary controls, i.e. , control is ef- 
fected through the data in a Neumann boundary condition. Let SI be a bounded open 
domain in ]R d , d = 2 or 3, and let T be its boundary. A simplified Ginzburg-Landau 
model for superconductivity is given by 

— At/h + (tpl + ipl + | A | 2 - l) Vh - V • (At/h) - A • Vt/> 2 = 0 in £2 , 

-At l) 2 + (ipl + t/>f + | A | 2 - l) t/> 2 + V • (At/q) + A • Vt/q =0 in £2 , 
n • (Vt/q + At/h) = ^9i onT, 

and 

n • (Vt/> 2 - At/q) = A# 2 on T . 

Here, Vh and Vh denote the real and imaginary parts, respectively, of the complex- 
valued order parameter, A is a given real magnetic potential, g\ and 32 are related 
to the normal component of the current at the boundary, and A > 0 is a “current 
loading” parameter. These equations are a special case of a more general model for 
superconductivity wherein A is also unknown; see, e.g., [22] for a derivation of the 
general model. It can be shown that in certain limits, e.g, high values of the applied 
Held, the above simpler model is valid; see [7]. 

By introducing appropriate rescalings, i.e, by replacing ipj by A ipj and gj by A gj, 
j = 1, 2, we can rewrite the above Ginzburg-Landau equations as follows: 


(4.25) 

-A Vh 

+ (|A| 2 - 

-l)^i 

- V • 

(At/ 2 ) - A • Vt / 2 + A (t / 2 

+ V 1 2 ) 

Vi = 0 

in £2 , 

(4.26) 

— At/h 

+ (|A| 2 - 

- 1)V’2 

+ v- 

(At/q) + A • Vt/q + A (t /> 2 

+ ^) 

to 

II 

O 

in £2 , 

(4.27) 



n • 

(VtG 

+ Aip 2 ) = A 31 on T , 




and 









(4.28) 



n • 

<] 

to 

- At/q) = A(/ 2 on T . 





We introduce the bilinear forms 

a(ip,<f>) = J ^Vt/> • V</> + ( | A | 2 — l)ip(pj dQ V ip,<p e H\n) 

and 

b(ip, <f>) = f A • (ip'Vtp — <pVip) d£l V ip, <f> E i7 1 (S2) . 

J a 

We assume that A £ H 1 (S2). Note that 

a(ip, <f) = a(<p, ip) and b(ip, <f>) = — b(<p, ip) . 
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Then, a weak formulation of the Ginzburg-Landau equations (4.25)-(4.28) is defined 
as follows: seek rf> = (t/q, ip 2 ) £ H 1 (S2) such that 

(4.29) (j>i) + b(ip 2 , <Pi) + A((V>i + ipl)ipi, </i) = Mgi, <Pi)r V ()i £ ff'fSl) 

and 

(4.30) a(ip 2 , fa) - Hipi, fa) + A ((-0i + i>l)fa, fa) = A(fif 2 , </ 2 }r V £ H 1 (ty . 

It can be shown that, for each g = (<7i , < 72 ) £ H _1 / 2 (r), (4.29) and (4.30) possess at 
least one solution ip G H 1 (S2) and that all solutions of (4.29) and (4.30) satisfy the a 
priori estimate 

(4-31) HV’illi + 1 1 '02 1 1 1 < C (| Iff 1 1 1 _ i/ 2 ,r + 1 1^2 1 1 — i/ 2 ,r) ; 

see, e.g., [11], for details. In the sequel, a solution of (4.25)-(4.28) will be understood 
in the sense of (4.29)-(4.30). 

Given a desired state ip 0 = ("0io , V’ 20 ) £ L 2 (I2), we define for any ip = (ip \ , fa) £ 
H 1 ^) and g = (< 71 , 172 ) £ L 2 (T) the functional 

(4.32) J(ip, g )=\ j ((Vh - V’io) 2 + (fa - fao) 2 ^) dn+^ j (g\ + g? 2 ) dT . 

We then consider the following optimal control problem associated with the Ginzburg- 
Landau equations for superconductivity: 

(4.33) min { J(ip, g) | ip £ H 1 (fl) , g £ 0 } subject to (4.29) and (4.30) , 
where 0 is a subset of L 2 (T). 

We define the spaces A = H 1 (S2), Y = (H 1 (S2))*, G = L 2 (T), and Z = 
[H 1 / 2+e (r2)]* where e £ (0, 1/2) is chosen such that H 1 (S2) H 1/,2+e (ff) L 4 (S2). 

By compact imbedding results, L 4 / 3 (S2) Z Y . For the time being, we assume 

that the admissible set 0 for the control g is a closed convex subset of G = L 2 (T). 

Let the continuous linear operator T £ £(Y;A) be defined as follows: for each 
f = (fijf'j) £ Y = (H 1 (S2))*, Tf = ip £ X = H 1 (S2) is the unique solution of 

a(fa,fa) + Kfa,fa) = (/1, V fa £ H\n) 


and 

a(fa, ^ 2 ) - Ki’i, fa) = (f 2 ,fa) v fa £ H 1 ^) . 

It can be easily verified that T is self-adjoint. Also, it can be shown that for most 
choices of A, the operator T is well defined; see [11]. 

We define the (differentiable) nonlinear mapping N : X —>■ Y by 


N(fa = 


(fal + fa)fa 
(fal + fa)fa 


V £ a 


or equivalently 

{N (ll>) , <j>) = ((-01 + 1pl)fa,fa) + ((fal + 1pl)fa,fa) V <t> = (fa, fa) £ A 
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and define K : H 1/,2 (T) — >■ Y as the injection mapping: 

(Kz, v) = — (z, v) r Vz £ H“ 1 / 2 (r),Vv £ H 4 (fi). 

Clearly, the constraint equations (4.29)-(4.30) can be expressed as 

if) + ATN(t/)) + XT Kg = 0 , 

i.e., in the form (2.2). With the obvious definitions for A(-) and £(•), i.e. , 

Hi>) = \J ((^ - ^10 ) 2 + (^2 - <M 2 ) dn 

and 

£{&)=\ Vg EG, 

the functional (4.32) can be expressed as 

JW,g) = a ^) + a 5(g), 

i.e., in the form (2.1). Thus, the minimization problem (4.33) is in the form of the 
minimization problem (2.3). 

We are now in a position to verify, for the minimization problem (4.33), all the 
hypotheses of §2 and §3. 

4-2.1. Verification of the hypotheses for the existence of optimal solutions. 

We first verify that the hypotheses (H1)-(H6) hold in the current setting. 

(HI) is obviously satisfied with a lower bound 0. 

(H2) holds with a = 1 and j3 = 2. 

(H3) is verified since if) = 0 and g = 0 is obviously a solution of (4.29)-(4.30). 

In order to verify (H4), we assume {g*-"- 1 } C 0 C L 2 (T) is a sequence satisfying 
g( n ) g i n L 2 (T); then, we have g( n ) g in H -1 / 2 (r) so that lim n ^oo (g*-"- 1 , v)r = 
(g, v )r for all v £ H 1 (S2), i.e., A'g( n ) A'g in Y. Assume that the sequence 
C H 1 (S2) satisfies if) in H 1 (S2); then, by using the compact imbed- 
ding H 1 ^) L 4 (fi), V> (n) V in L 4 (fi). Now, 

+ ( 4 n ) ) 2 )^ n ) ^ i ) + (((^ n) ) 2 + 2 ) 

-► +^ 2 )^ 2 , ^ 2 ) = {N(rl>),<j>) . 

Hence, (H4) is verihed. 

The verification of (H5) follows directly from the observation that the mappings 
4 ^ F(4) = (l/2)||^ - V>o||o and g >->■ 5(g) = (l/2)||g|| 2 r are convex. 

The verification of (H6) is a trivial consequence of the a priori estimate (4.31). 

It is now just a matter of citing Theorem 2.1 to prove the existence of an optimal 
solution that minimizes (4.32) subject to (4.29)-(4.30). 

THEOREM 4.4. There exists a (<j>, g) £ H 1 (S2) x 0 such that (4.32) is minimized subject 
to (4.29)-(4.30). □ 

4-2.2. Verification of the hypotheses for the existence of Lagrange multipliers. 
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We now assume (ip, g) is an optimal solution and turn to the verification of 
hypotheses (H7)-(H9). 

The validity of (H7) is obvious. 

(H8) holds since the mapping g h- > £(g) = (1/2) J r |g| 2 dr is convex. 

(H9) can be verified as follows. For any ip E X , the operator N'(ip) : X —>■ Y is 
given by 


N'(ip)-<f> 


(3-ipl + i + (2ip 1 ip 2 )<p2 
(3i Pl + ipl)<p 2 + (2ipiip2)<pi 


V <p — (<f> i, (p 2) E X . 


Thus, we obtain that N'(ip) • <j> E L 4 / 3 (fi) ^ [H 1 / 2+£ (fi)]* = Z. 
The Lagrangian is given by 


=kj(ip, g) 

- m) + Ki’ 2 , + a((v>i + ipl)ipi, m) - Mgi,m)r 

+ a(ip 2 ,r] 2 ) - - A ((/ , 1 + i’ 2 )i’ 2 ,V 2 ) - A( 02 , j? 2 )r} 

for all (0 g,i),l:)eIxGxFxR = H 1 (S2) x L 2 (T) x H 1 (S2) x R . Note that in this 
form of the Lagrangian, the Lagrange multiplier r) E X = Y* so that we have already 
introduced the change of variables indicated between (2.17 )-(2 .18) and (2 . 19)-(2 .21). 

Having verified the hypotheses (H7)-(H9), we may apply Theorem 2.4 to conclude 
that there exists a Lagrange multiplier r) E X = H 1 (S2) and a real number k such that 

(4.34) r) + XT* ([N'(iP)}* -n-kJ^iP, g)) = 0 
and 

(4.35) C(ip,g,r),k) < C(ip,z,r),k) V z E 0 


and that for almost all values of A, we may choose k = 1. 

Recall that T is self-adjoint. Also, note that for any ip E X = H 1 (S2), 


, _ ( (3-0? + ipl)rn + (20i 02 ) ?72 

V (30i + ZOm + (20 i 0 2 )»?i 


[N'm* -v= rZa T 2z : (*,*)€*. 


Thus, N'(ip) is self-adjoint as well and (4.34), with k = 1, can be rewritten as 

-&(Ci,d2) + A((3t/>( + 0l)?7i, Cl) 


(4.36) 
and 

(4.37) 


+ a((20 1 0 2 )t 72 ,Ci) = A (01 - 010 , Cl) V Cl e H\n) 

a(C 2,m) + H( 2,m) + A ((30l + 0i02, C 2 ) 

+ A((20!0 2 )»7i, C 2 ) = M’4’2 — 020 , C 2 ) V C 2 E 77 1 (fl) . 


Using the definition of the Lagrangian functional, (4.35), with k = 1, can be rewritten 
as 

- (z, z)r + A(z, 0 r - 77 (g, g)r - A(g, 0r > 0 Vz£0. 
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Note that, in the above expression, we have already employed hypothesis (H12) which 
in the current context is trivially satisfied with E the identity operator on G* = G = 

L 2 (r). For each e £ (0, 1) and each t £ 0, set z = et + (1 — e)g £ 0 in the last 

equation to obtain 

e 2 

Y (* — S, * — g)r + e (t — g, g)r + e (t — g, rf) r > 0 V t £ 0 

so that, after dividing by e > 0 and then letting e —>■ 0 + , we obtain 

(4.38) (t - g,g + j?) r > 0 V t £ 0 . 

We see that for almost all values of A, necessary conditions for an optimum are 
that (4.29)-(4.30) and (4.36)-(4.38) are satisfied. Again, the system formed by these 
equations will be called an optimality system. 

We now specialize to the case 0 = L 2 (T). Note that the hypothesis (H10) is 
satisfied. Then, using Theorem 2.5, we see that the inequality (4.38) becomes an 
equality and, by letting z = t — g vary arbitrarily in L 2 (T), we now have, instead of 

(4.38) , 

(4.39) (z, g + q) r = 0 Vz£L 2 (r). 

Thus, according to that theorem, we have that for almost all A, an optimality system 
of equations is now given by (4.29)-(4.30), (4.36)-(4.37), and (4.39). However, we can 
go further and verify that the hypothesis (Hll) is valid, which in turn will justify the 
existence of a Lagrange multiplier satisfying the optimality system for all A £ A. 

To verify (Hll), we first note that, through the change of variable £ = T*v, that 
assumption can be equivalently stated as follows: 

if £ £ Y* satisfies (i + XT* [^'(m)]*)^ = 0 and K*£ = 0, then £ = 0 . 

To verify this version of (Hll), we assume that (, ) ' H 1 (S2) satisfies (/ + 

XT* [A r '(V’)]*)£ = 0 and K*£ = 0, i.e., 

a(Ci,Ci) — KCi) £ 2 ) + ^ ((3i/ , i + t/’D'fi) Ci) 

+a((2^ 2 )6,Ci) = 0 V(i en 1 ^), 

a (( 2 , £ 2 ) + b(( 2 , Cl) + ^ ((3C , 2 + "01 )<? 2 , C 2 ) 

+A ((20! 0 2 )<f r , C 2 ) = 0 VC 2 ei0(«), 

and 

^ = 0 on r. 

(Note that K*£ = C|r-) Let fi be a smooth extension of £2 such that £2 is a compact 
subset of £2 . We then define (, , if) and A' to be the extension, by zero outside £2, 
of if) and A, respectively. Let the forms a'(-, •), h (•, •), and (•, •)' defined over £2 
be the analogues of corresponding forms defined over £2. We may show from the last 
three equations that 

t' eH'fSl'), if>' £ L 6 (£2 < ) , 

a'(Ci,0) - b\( uQ + X ((3^ 2 + 0> 2 )CUi)' 

+A((2^)£,Ci)' = 0 VCieHj(£2'), 
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and 


a '(C'2, £2) + ^ (C2, 1 ) + ^ (( 3 V 4 + V 4 ) < s2 > C2) 

+A((2^^)^,C 2 )' 


In the sense of distribution, satisfies 


0 V( 2 eHj(Q'). 


(4.40) 


A0 - 2A' • V0 + (|Af + A (3-01 2 + ^ 2 ) - 1)£ 

-(V • A' - 2A 0102)^2 = 0 in ft' 


and 

(4.41) 


A& + 2A' • V0 + (V • A' + 2A 00' )£' 

+ (|A'| 2 + A (302 + i’l ) ~ l )^2 = 0 in S2 . 


We now quote the following unique continuation result whose proof can be found in 
[17]. See also [12] and [19]. 

LEMMA 4.5. Let L2 be an open and connected subset of ]R d , d = 2 or 3. Let the 
functions V £ 00(12 )] dxd for some q > 2 and W £ [L 1 2 0 d c _1 (S2 )] dxdxd g tve n. Lf 
l e H00') , -A fi + Y?j=i J2k=i Wijkidfk/dxj) + Y?j=i Vijtj = 0 (in the sense of 

distributions), i = 1, and £ = 0 on an open, non-empty subset of Ll , then £ = 0 

on LI . □ 

Since A £ H 1 (S2) and rf> £ H 1 (S2), it is easy to see that the coefficients in (4.40)- 
(4.41) satisfy the regularity requirements of Lemma 4.5. Also note that ^ = 0 on 
(LI \ 12) which contains an open set. Thus we obtain that (, = 0 in 12 , or ^ = 0 in 12, 
i.e. , (Hll) is verified. 

Hence we conclude that for all A, the optimality system (4.29)-(4.30), (4.36)- 
(4.37), and (4.39) has a solution. Thus, we have Theorem 2.6 which, in the present 
context, is given as follows. 

THEOREM 4.6. Let (xj), g) £ H 1 (S2) x L 2 (r) denote an optimal solution that minimizes 
(4.32) subject to (4.29)-(4.30). Then, for all A £ A, there exists a nonzero Lagrange 
multiplier r) £ H 1 (S2) satisfying the Euler equations (4.36)-(4.37) and (4.39). □ 


J t .2.3. Verification of the hypotheses for approximations and error estimates. 

We finally verify the hypotheses (H13)-(H19) that are used in connection with 
approximations and to derive error estimates. 

A Unite element discretization of the optimality system (4.29)-(4.30), (4.36)- 
(4.37), and (4.39) is defined in the usual manner. We first choose families of Unite 
dimensional subspaces X h C 10(12) and G h C L 2 (T) parameterized by a parameter h 
that tends to zero and satisfying the following approximation properties: there exists 
a constant C and an integer r such that 


(4.42) inf ||^-^||i < Ch m \\4>\\ m+ i , V* £ H m+1 (12), 1 < m < r 

<f> h £X h 


and 


(4.43) 


inf llz — z ft ||or < Ch m 

z h £G h ' vgH 


inf 

+ 1 / 2 (n), v| r =z 


IML + 1/2 , 


Vz £ H m+1/2 (12)| r , 1 < m < r. 
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One may consult, e.g., [8] and [15], for some finite element spaces satisfying (4.42) and 

(4.43) . For example, one may choose X h = V h x V h where V h is the piecewise linear 
or quadratic finite element space defined with respect to a family of triangulations 
of £2. In this case, h is a measure of the grid size. For simplicity we may choose 
G h = (X ft )|r, i.e. , the functions in G h are the restrictions to the boundary T of 
functions belonging to X h . 

Once the approximating spaces have been chosen, we may formulate the approx- 
imate problem for the optimality system (4.29)-(4.30), (4.36)-(4.37), and (4.39): seek 
ij> h £ X h , g h £ G h , and T) h £ X h such that 

(4.44) a(0, <) + b« <) + A{ ((^) 2 + OYWiXi) = X(gl <)r V < £ V h , 

(4.45) a(il >2 , < 2 ) - <4, 4) + M (Of + OfO , 4 ) = <92,4 r V < 2 £ V h , 

<Ci, 4) - biCtri) + A((3 Of + OfOiX 1 ) 

+ A ((2^^)^, Ci) = <4 ~ rhoXi) V < G V h , 

<(2,4) + <44) + HAOf + Of)A,4 

+ x((244<,4 = «2 - <o,<) v (2 e , 


(z\g ft +J? ft )r = 0 Vz h £G h . 

The operator T h £ C(Y;X h ) is dehned as follows: for f £ Y, T h f = tf) h £ X h is 
the solution for 

aO,4) + «2,4) = (h,4) V<£V h 

and 

a(0,<2)-« h i,4) = {f 2 ,< 2 ) v4ev h . 

Since T = T* , we define ( T*) h = T h . 

We dehne the operator E h : L 2 (T) — >■ G h as the L 2 (r)-projection on G h , i.e., for 
each g G L 2 (T), 

(E h g,z ft )r = (g,z ft ) r Vz h £ G h . 

Since G = L 2 (T) is reflexive, E h is in fact an operator from G* — >■ G h . 

By results concerning the approximation of the Ginzburg-Landau equations (see, 
e.g., [11]), we obtain 

||(T — T h )f\\x —>■ 0 

as h —>■ 0, for all f G Y . This is simply a restatement of (H13). 

(H14) follows trivially from (H13) and the fact that T is self-adjoint and we have 
chosen (T*) h =T h . 

(H15) follows from the best approximation property of L 2 (r)-projections and 
(4.43). 

(H16) and (H17) follow from the fact that N and T are polynomials. Here we 
also use imbedding theorems and Cauchy inequalities. 


(4.46) 


(4.47) 


(4.48) 
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Setting Z = Z = H 1/,2+e (fl), we have that Z [H 1 (S1)]* = X* . For each 

H 1 ^) and £ G H 1 (S1), Sobolev imbedding theorems imply that 




mi + m)m + mm)m 
mi + + mm)m 


E L 4 / 3 (Q) cl, 


([^'W •<)••? 


(6"0lCl + 2V> 2 C2)»7l + mmm + (2Cl-02)??2 \ -,-4/3/ro r 7 

( 6 ^ 2(2 + 2V>iCi)»72 + (2V>iC2)»7i + (2CiV , 2)»7i / ^ 


and 

mmo v= e l 4 / 3 (q) c z . 

These relations verify (H18). 

From the definition ofthe operator K we see that K maps L 2 (T) into [H 1/,2+e (£!)]* , 
i.e. , K maps G into Z. Thus (H19) is verified. 

Hence, we are now in a position to apply Theorem 3.5 to derive error estimates 
for the approximate solutions of the optimality system (4.29)-(4.30), (4.36)-(4.37), and 
(4.39). It should be noted that Lemma 3.4 implies that for almost all values of A, the 
solutions of the optimality system are regular. 

THEOREM 4.7. Assume that A is a compact interval of 1R_|_ and that there exists a 
branch {(A, ■0( A) , g(A), 77(A)) : A £ A) of regular solutions ofthe optimality system 
(4.29)-(4.30), (4.36)-(4.37), and (4.39). Assume that the finite element spaces X h and 
G h satisfy the hypotheses (4.42)-(4.43). Then, there exists a 6 > 0 and a ho > 0 such 
that for h < ho, the discrete optimality system (4.44)-(4.48) has a unique branch of 
solutions {(A, xj} h (\), g h (X), rj h (X)) : A G A} satisfying 


{||^(A) - V-(A)||i + ||g ft (A) - g(A) 1 1 o,r + ||»? ft (A) - »?( A) 1 1 x 5 for all A G A. 


Moreover, 

lim{||^(A) - V>(A)||! + ||g ft (A) - g(A)|| 0 r + ||^(A) - »?( A) 1 1 x } = 0 , 

^ 0 


uniformly in A £ A. 

If, in addition, the solution ofthe optimality system satisfies (^>(A), g(A), jj(A)) G 
H m+1 (S2) x H m+1 / 2 (S2)|r x H m+1 (S2) for A G A, then there exists a constant C, 
independent of h, such that 

IIV’(A) - t A) 1 1 i + ||g(A) - g ft (A) | | 0 ,r + ||i?(A) - ^(A)|U 

<Ch m ( ||V’(A)||m + l + inf || v ||m + l/2 + ||*?(A)||m + l) , 

vGH”*+ 1 /2(n),v| r =g 

uniformly m A £ A. 

Proof: All results follow from Theorem 3.5. For the last result, we also use (3.25) and 
the estimates (see [11]) 

|| mm - 7)V>||i < Ch m \m m+1 for E H m+1 (fl) , 

||((T*) h (T*)-i - I)n 111 = IKT^T - 1 - mw, < C/i m ||jj|| m+1 for r, G H^fi) , 
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and (see, e.g., [2], [8] and [15]) 


\\{E h E - 1 


/)g||o,r < Ch n 


inf ||v|| m+ i/ 2 

veH”»+ 1 / 2 (n),v|r=g 


for g E H m+1 / 2 (£!)| r • 


In these estimates, the constant C is independent of h, if), g, jj, and A. □ 
Remark. In fact, we obtain from (4.39) that g = — *?|r which implies 


inf 

v GH "*+ 1 / 2 (n),v|r=; 


|| v ||m + l /2 < ||*?||m + l /2 ^ ll*?l|m + l 


so that the term (inf vgH m+i/ 2 ( n ) v | r =g ll v llm+i/ 2 ) in the right-hand side of the error 
estimate is redundant. □ 

Remark. By using (4.39) again, along with (4.48) and the error estimate in Theorem 
4.7, we have the following improved error estimate for the approximation of the control 

g: 


||g(A) - g ft (A)|| 1/2ir < C\\rj(X) - ^(A)^ < Ch m ( ||^(A)|| m+1 + ||»?(A)|| m+1 ) . 

Of course, we also use the fact that we have chosen G h = (X ft )|r C H 1 I 2 (Y). □ 

4.3. Dirichlet boundary control for the Navier-Stokes equations of incom- 
pressible, viscous flow 


For this application we will use Dirichlet boundary controls, i.e. , control is effected 
through the data in a Dirichlet boundary condition. Let SI denote a bounded domain 
in IR d , d = 2 or 3 with a boundary denoted by T. Let u and p denote the velocity and 
pressure Helds in SI. The Navier-Stokes equations for a viscous, incompressible flow 
are given by (see, e.g., [13], [14], or [20]) 


-v V • ((Vu) + (Vu) T ) + (u • V)u + Vp = f in SI , 


V • u = 0 in £2 


and 

u = b + g on T, 

where f is a given body force, b and g are boundary velocity data with f r b • n dr = 0 
and J r g • n dr = 0, and v denotes the (constant) kinematic viscosity. We have 
absorbed the constant density into the pressure and the body force. If the variables in 
these equations are nondimensionalized, then v is simply the inverse of the Reynolds 
number Re. 

Setting A = 1/v = Re and replacing p with p/ A, b with Ab, and g with Ag, we 
may write the Navier-Stokes equations in the form 

(4.49) -V • ((Vu) + (Vu) T ) + Vp + Au • Vu = Af in SI , 


(4.50) V • u = 0 in SI, 
and 

(4.51) u = A(b + g) on f. 
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We introduce the subspaces 


L 2 0 (n) = 


p E L 2 (Cl) 


p di2 = 0 


and 


H „(r) = jg e H^r) 


I g 

r 



We also introduce the bilinear forms 

a(u, v) = t / ((Vu) + (Vu) T ) : ((Vv) + (Vv) T ) dQ V u, v E H^fi) 
^ J n 

and 

&(v, q) = — qV ■ v di2 VvE H 1 (S1) and Vp£ L 2 (i2) 

J a 

and the trilinear form 


c(u,v,w)= / (u-V)v-wdSl V u, v, w E H 1 (S1) . 

Jn 

These forms are continuous over the spaces of definition indicated above. Moreover, 
we have the coercivity properties 


(4.52) 
and 

(4.53) 


a( v , v) > C a ||v||? V v £ H((Q) 


sup > Olkllo V g G i„(Q) 

o^veHj(n) ll v lli 


for some constants C a and Cj > 0. For details concerning the notation employed 
and/or for (4.52)- (4.53) , one may consult [13], [14], and [20]. 

We recast the Navier-Stokes equations (4.49)-(4.51) into the following particular 
weak form (see, e.g., [15]): seek (u,p>, t) E H 1 (S2) x T 2 (£2) x H -1 / 2 (r) such that 

(4.54) a(u, v) + J(v, p) - (t, v) r + A c(u, u, v) = A (f , v) VvEH'ffi), 


(4.55) b(u,q) = 0 V?Gi^(0), 
and 

(4.56) (s,u) r - A (s,g) r = A (s,b) r V s E H _ 1 / 2 (r) . 

Formally we have 

t = [-pn + (Vu + (Vu) T ) • n] r , 

i.e. , t is the stress force on the boundary. The existence of a solution (u,p, t) for the 
system (4.54)-(4.56) was established in [15]. 
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Given a desired velocity field uo, we define for any (u,p, t) G H 1 (S1) x L 2 (£l) x 
H-i/ 2 (r) and g G H^(T) the functional 

(4.57) J(u,p,t,g) = j jju-u 0 \ 4 dn+^ j^ (\V s g\ 2 + \g\ 2 )dT , 
where V s denotes the surface gradient. 

We define the spaces X = H 1 ^) x L 2 0 { Q) x H” 1 / 2 ^), Y = [H 1 ^)]* x Lg(ft) x 
H 1 / 2 (r), G = Hjj(r), and Z = L 3 / 2 (S1) x {0} x H 1 (r). By compact imbedding results, 
Z is compactly imbedded into Y . For the time being, we assume that the admissible 
set 0 for the control g is a closed, convex subset of G = H^T). 

We then consider the following optimal control problem associated with the 
Navier-Stokes equations: 

(4.58) min{ J(u, p, t , g) : (u, p, t) 6 1, g 6 0} subject to (4.54)-(4.56) . 

We define the continuous linear operator T G £(Y;X) as follows: for each 
G Y, T(C,t),k ) = (ri , p, t) G X is the unique solution of 

a(u, v) + b(v,p) - (t, v) r = ((, v) V v G H^fi) , 

b(u,q) = V ? £ f^(Q) 

and 

(s,u) r = (s,/c) r V s G H _1/2 (r) . 

It can be easily verified that T is self-adjoint. 

We define the (differentiable) nonlinear mapping N : X —>■ Y by 

/f-u Vu 

N(u,p,t) = - 0 

V b 

or equivalently 

{N(u,p,t), (v,q, s)) = -(f, v) + c(u,u, v) - (s,b) r V (v,q, s) G X 
and define K : H 1/,2 (r) Y by 


Kg 


or equivalently 

(Kg, (v,q, s)) = — (s, g) r V g G H 1 / 2 (r) ,V (v,q, s) G X . 
Clearly, the constraint equations (4.54)-(4.56) can be expressed as 
(u ,p, t) + \ TN(u,p, t) + XT Kg = 0 , 
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i.e., in the form (2.2). With the obvious definitions for T{-) and £(•), i.e. , 


J r (u,p,t)=j f |u — u 0 1 4 dQ V(u,p,t)Gl 
4 Jn 

and 

£(g) = \ / r (|v sg | 2 + |g| 2 ) dr, 

the functional (4.57) can be expressed as 


J(u,p,t,g) = AX(u,p,t) + A 5(g) , 


i.e., in the form (2.3). 

We are now in a position to verify, for the minimization problem (4.58), all the 
hypotheses of §2 and §3. 

4-3.1. Verification of the hypotheses for the existence of optimal solutions. 

We first verify the that the hypotheses (H1)-(H6) hold in the current setting. 

(HI) is obviously satisfied with a lower bound 0. 

(H2) holds with a = 1 and /? = 2. 

(H3) is verified with the choice (u^ 0 / p(°\ t(°), 0) £ X x 0 where (u( 0 /p( 0 )) is 
a solution to the Navier-Stokes equations with Dirichlet boundary conditions, and 
t(°) = [-p(°) n + (Vu<°) + (Vu(°)) T ) • n] r ; see, e.g., [13] or [20], 

In order to verify (H4), we assume {g*-"- 1 } C 0 C H/(T) is a sequence satisfying 
g( n ) g i n H 1 (r); then we have g( n ) g in H 1/,2 (T) so that lim n ^oo(g( n ), v)r = 
(g, v)r for all v £ H 1 (S2), i.e., A'g( n ) Kg in Y . Assume that the sequence {u( n )| C 
H 4 (ft) satishes u( n ) u in H 1 (S2); then u( n ) — >■ u in L 4 (S2) by the compactness of 
the imbedding H 1 (S2) L 4 (S2). Now, 

(N( u^), v) = c(u( n \ u( n \ v) = c(u, u( n \ v) + c(u( n ) - u, u( n \ v) 

->■ c(u, u, v) + 0 = (N( u) , v) as n —>■ 00 . 


Hence, (H4) is verihed. 

The verihcation of (H5) follows directly from the observation that the mappings 
(u,p,t) T(u,p,t) = (1/4) ||u — u||L4(n) and g ^ £(g) = (1/2) ||g||yr are convex. 

To verify (H6), we combine a priori estimates obtained from the constraint equa- 
tions and the functional. Let p( k \ t^ k \ g*-*- 1 } C H 1 (S2) x T 2 (f2) x H _1 / 2 (r) x 

H /(r) be a sequence such that 

(4.59) J(u (fc) ,g (fc) ) < C, 

(4.60) afuW.vj + Hv^^-^/vJr + Acfu^/uWv) = A (f , v) Vv e H 1 (fi), 

(4.61) b(u( k \q) = 0 V?ei](0), 
and 

(4.62) (s,u( fc )) r - A(s,g( fc )) r = A(s,b) r Vser^f). 
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First, (4.59) implies that (u( fc ),g( fc )) is uniformly bounded in L 4 (S1) x H 1 (T). For each 
we may choose a £ H 1 (S1) x Tq(£ 1) that satisfies the Stokes problem 

(4.63) a(w (k \v) + b(v,r (k y) = (i,v) Vv£H)(fi), 


(4.64) b(w (k \q) = 0 Vg£Lg(fi) 
and 

(4.65) w*-^ = A(g ( - fc - ) + b) on T. 

Furthermore, there holds the estimate 

(4-66) l| W ( i) l| 1 <C'(||f||o+||b|| 1/2 , r + 11^11^). 

By subtracting (4.63) from (4.60) with v = — w( k \ also using (4.61) and (4.64), 

we obtain 


(4.67) 


a(u ( - fc - ) — w( k \ u*-^ — w*-^) 


-AcfuW^W^W-w^) 

AcfuW.uW-w^^W). 


Note that 


|c(u^,u^- w^,u^)| 
_ 1 
~~ 2 

< C 

1 


f u (k ~> ■ ((V(u^ - w (k '>)) + (V(u^ - w^)) T ) • u (k ~> <m 

J a 


< 


4A 


(V(u^ - w^)) + (V(u^ - w^)) T 

(V(u^ - w^)) + (V(u^ - w ( fc ))) T 



u(*) 

0 



■C x 


L 4 (n) 


L 4 (n) 


so that, using (4.67), we have that 


(V(u^ - w (k '>)) + (V(u^ - 


< C A 


,(*) 


L 4 (n) 


Then, by (4.66) and the triangle inequality, we have that 


(Vu (i) ) + (Vu<*)) T < C{||f ||o + ||b|| 1/ 2 ,r + ||g (fc) ||i,r + l|u^||^ (n) } . 


Thus, 


||(Vu( fc )) + (Vu( fc )) T ||o + ||u( fc )|| 0 ,r 
< ||(Vu( fc )) + (Vu( fc )) T ||o + ||b|| 0 ,r + ||g (fc) ||o,r 
<C7(||f||o + ||b|| 1/2ir + ||gW|| lir +||uW||2 4(n) ). 

Since the mapping u h- > ||Vu + (Vu) T ||o + ||u||o,r defines a norm on H 1 (S1) equivalent 
to the standard H 1 (f2)-norm, we have that 

l|u (fc) ||i < C{||f ||o + ||b|| 1/ 2 ,r + ||g (fc) ||i,r + l|u^||L (n) } , 
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and, since Hu^Hl^) and ||g ( - fc ' ) ||i 1 r are uniformly bounded, we conclude that ||u ( - fc )||i 
is uniformly bounded as well. One easily concludes from (4.60) that | |tt & ) 1 1 _ x / 2 ,r is 
uniformly bounded. Thus (H6) is verified. 

It is now just a matter of citing Theorem 2.1 to conclude the existence of an 
optimal solution that minimizes (4.57) subject (4.54) — (4.56). 

THEOREM 4.8. There exists a (u,p, t,g)E H 1 (S2) x L 2 (Ll) x H -1 / 2 (S2) x O such that 
(4.57) is minimized subject to (4.54)-(4.56). □ 

4-3.2. Verification of the hypotheses for the existence of Lagrange multipliers. 

We now assume (u,p, t,g) is an optimal solution and turn to the verification of 
hypotheses (H7)-(H9). 

The validity of (H7) is obvious. 

(H8) holds since the mapping z h- > £(g) = (1/2) / r (|V s g| 2 + |g| 2 ) dr is convex. 
(H9) can be verified as follows. For any (u,p, t) E X, the operator N'(u,p, t) : 
X — >■ Y is given by 


N'(u,p,t) ■ (v,q, s) 


u • Vv + v • Vu 
0 
0 


for all (v,q, s) E H 1 ^) x L 2 0 (fl) x H“ 1 / 2 (r). Thus we obtain N'( u, p, t) • (v, q, s) E 
L 3 / 2 (S2) x {0} x H^r) = Z. 

The Lagrangian is given by 


£(u,p,t,g,i/, <f>,T, k) 

= kj( u, g) - {a( u, v) + A c(u, u, v) + b{v, p) + b( u, <j>) - (r, u) r 

- (t, v) Y ' - A (f , v) Y ' + A (r, b) r + A (r, g) r } 

for all (u, p, t, g, v, <f>, t, k) eX xGxX x 1R = H^fi) x L 2 0 (Q) x H“ 1 / 2 (r) x H/(r) x 
H 1 ^) x L 2 (Ll) x H“ 1 / 2 (r) x 1R. Note that in this form of the Lagrangian, the 
Lagrange multiplier (v, <j>, r) E X = Y* so that we have already introduced the change 
of variables indicated between (2.17)-(2.18) and (2 . 19)-(2 .2 1) . 

Having verified the hypotheses (H7)-(H9), we may apply Theorem 2.4 to conclude 
that there exists a Lagrange multiplier ( v,<j>,T ) E X = H 1 (S2) x L 2 (^) x H -1 / 2 (r) 
and a real number k such that 

(4.68) {v,<f>,r) + XT* ([N'(u,p,t)]* ■ {v,<f>,r) - k J'( UiPit )(u, p, t, g)) = 0 
and 

(4.69) C(u,p,t,z,v,<f>,T,k) < C(u,p,t,g,v,<f>,T,k) VzE0 

and that for almost all values of A, we may choose k = 1. 

Recall that T* = T . Also, note that for (u, p, t) E X = H 1 (S2) xL 2 (Ll) xH -1 / 2 (r), 
the operator [tV'(u, p, t)]* : X —>■ Y is given by 

/ -u • Vv + v • (Vu) T \ 

[tV'( u ,p,t)]* • (v, ? ,s) = f 0 1 V (v, g, s) £ X. 
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Thus, (4.68), with k = 1, can be rewritten as 


a ( w ; v) + \ c(w, u, v) + A c(u, w, v) + b{ w, <^) - (r, w) r 

= A ((u - u 0 ) 3 , w) VwEH 1 (fi), 

(4.71) b(v, r) = 0 VreTg(fi), 
and 

(4.72) (y,v)r = 0 VyeH- 1 ^). 

In the right-hand side of (4.70), we use the notation (v 3 , w) = J^_i(i> 3 , wj). 

Using the definition of the Lagrangian functional, (4.69), with k = 1, can be 
rewritten as 


^ (V s z, V s z) r + ^ (z, z) r - ^ (V s g, V s g) r 

-■^■(g)g)r- A(r,z)r + A(r,g) r > 0 VzE0. 

For each e E (0, 1) and each z E 0, by plugging ez + (1 — e)g E 0 into the last 
inequality we obtain 

e 2 

e (V.g, V s (z - g)) r + e (g, z - g)r + y (V s (z - g), V s (z - g)) r 

e 2 

+ Y ( z “ g’ z “ g)r - e ( r > z - g)r > 0 V z E 0 

so that, after dividing by e > 0 and then letting e —>■ 0 + , we obtain 

(4.73) (V s g,V s (z-g)) r + (g,z-g) r - (r,z) r > 0 VzE0. 

We see that for almost all values of A, necessary conditions for an optimum are 
that (4.54)-(4.56), (4.70)-(4.72) and (4.73) are satisfied. Again, the system formed by 
these equations will be called an optimality system. 

We now specialize to the case 0 = H^(r). Note that the hypothesis (H10) is 
satisfied. Then using Theorem 2.5, we see that the inequality (4.73) becomes an 
equality and, by letting z = k — g vary arbitrarily in H*(r), we now have, instead of 

(4.73) , 

(4.74) (V s g,V s z)r + (g,z) r -(r,z)r = 0 Vz 6 0 = Hi(r). 

Thus, according to that theorem, we have that for almost all A, an optimality system 
of equations is now given by (4.54)-(4.56), (4.70)-(4.72) and (4.74). However, we can 
go further and verify that the hypothesis (Hll) is valid, which in turn will justify the 
existence of a Lagrange multiplier satisfying the optimality system for all A £ A. 

We now verify (Hll) which we again note can be equivalently stated as follows: 

if £ E Y* satisfies (i + XT* [77'(m)]*)^ = 0 and K*f = 0, then £ = 0 . 
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To verify this hypothesis, we assume that (£, a,0) E Y* = H 1 (S1) x T 2 (£2) x H 1/,2 (T) 
satishes (/ + A T* [At'(u, p, t)]*) (£, a, 0) = (0,0, 0) and K*(£, cr,0) = 0, i.e. , 


a(w,£) + A c(w, u,£) + A c(u, w,£) + b( w, u) - ( 6 , w) r = 0 V w £ H 1 (fi) , 

&(C,r) = 0 Vrei^), 

(y,C)r = o VyeH-^r), 


and 

0 = 0 on r. 

(Note that K*(£,a,0) = 0 .) Let SI be a smooth extension of SI such that SI is a 
compact subset of SI . We then set (, , a’ and u' to be the extension, by zero outside 
SI, of a and u, respectively. We may show from the last four equations that 

{'eH 1 (fi'), <r'eLl( n'), 

(4.75) a'(w,£') + A c'(w, u',£') + A c'(iT, w,£') + b (w, cr') = 0 V w E Hj(S2 ) , 
and 

(4.76) b\t',r) = 0 V r e L 2 ( Si'), 

where the forms a'(-, •), b (•, •) and c'(-, •, •) dehned over SI are the analogues of cor- 
responding forms dehned over SI. Using a unique continuation result for the system 
(4.75)-(4.76) that was established in [16] or [17], we obtain = 0 and cr' = 0 in SI , 
or ^ = 0 and a = 0 in SI. Thus (Hll) is verified. 

Hence we conclude that for all A, the optimality system (4.54)-(4.56), (4.70)- 
(4.72), and (4.74) has a solution. Thus, we have Theorem 2.6 which, in the present 
context, is given as follows. 

Theorem 4.9. Let (u,p,t,g) e H 1 (S)) x U 2 (S2) x H _1 / 2 (r) x H* (r) denote an optimal 
solution that minimizes (4.57) subject to (4.54)-(4.56). Then, for all A E A, there 
exists a nonzero Lagrange multiplier ( v,<j>,T ) E H 1 (S2) x T 2 (£2) x H -1 / 2 (r) satisfying 
the Euler eguations (4.70)-(4.72) and (4.74). □ 

Note that, in the above expression, we have already employed hypothesis (H12) 
which in the current context is easily seen to be satisfied with E : G —>■ G* dehned by 

(Eg, z) = ^(V s g • v s z + g • z) dr v z E Hi (r) = g . 

We also note that for each hxed r, (4.74) with g E H* (r) is equivalent to 

(4.77) (V s g,V s k)r + (g,k) r + 7 ^k-ndr = (r,k) r V k E H 1 (r) 

and 

(4.78) J g ndT = 0, 


40 



where 7 £ ]R is an additional unknown constant that accounts for the single integral 
constraint of (4.78). The equivalence can be shown as follows. First, an application 
of Lax-Milgram Lemma to (4.74) on the space H* (T) guarantees the existence and 
uniqueness of a solution g £ H*(r) to (4.74); this solution g clearly satisfies (4.77)- 

(4.78) with 7 = Jp (t • n — V s g : V s n — g • n) dr. Conversely, any solution (g,y) of 
(4.77)-(4.78) trivially satisfies (4.74). Although (4.74) and (4.77)-(4.78) are equivalent, 
the latter is more easily discretized. 

4-3.3. Verification of the hypotheses for approximations and error estimates. 

We finally verify the hypotheses (H13)-(H19) that are used in connection with 
approximations and error estimates. 

A finite element discretization of the optimality system (4.54)-(4.56), (4.70)- 
(4.72), and (4.74) is defined as follows. First, one chooses families of finite dimensional 
subspaces C H 1 (S2) and S h C L 2 (£2). These families are parameterized by the pa- 
rameter h that tends to zero; commonly, this parameter is chosen to be some measure 
of the grid size in a subdivision of £2 into finite elements. We let Sq = S h C L 2 (£2) and 
V([ = V h nHj(£2). 

One may choose any pair of subspaces and S h that can be used for find- 
ing finite element approximations of solutions of the Navier-Stokes equations. Thus, 
concerning these subspaces, we make the following standard assumptions which are 
exactly those employed in well-known finite element methods for the Navier-Stokes 
equations. First, we have the approximation properties: there exist an integer k and 
a constant C , independent of h, v and q, such that 

(4.79) inf ||v — v ft ||i < C/i m ||v|| m+ i V v £ H m+1 (£2) , 1 < m < k 

■yh 

and 


(4.80) inf \\q-q h \\o < Ch m \\q\\ m V q £ C Lg(£2) , 1 < m < k . 

q h es J 

Next, we assume the inf- sup condition, or Ladyzhenskaya-Babuska-Brezzi condition: 
there exists a constant C , independent of h , such that 


(4.81) 


inf 

0 ?q h €S*; 


sup 


b(v h ,q h ) 


>c. 


This condition assures the stability of finite element discretizations of the Navier- 
Stokes equations. For thorough discussions of the approximation properties (4.79)- 
(4.80), see, e.g., [2] or [8], and for like discussions of the stability condition (4.81), see, 
e.g., [13] or [14]. The latter references may also be consulted for a catalogue of finite 
element subspaces that meet the requirements of (4.79)-(4.81). 

Next, let = V ft |r, i.e. , consists of the restriction, to the boundary T, of 
functions belonging to ~V h . For all choices of conforming finite element spaces V ft , e.g., 
Lagrange type finite element spaces, we have that P ft C H - 1 / 2 (r). For the subspaces 
P ft = V ft |r, we can show the following approximation property: there exist an integer 
k and a constant C , independent of h and s, such that 


inf \\s-s h 


s H-1/2, r 


< Ch m inf 1 1 v I 

vEH m (n),v|r=s 


V s £ H m (£2)| r , 1 < m < k. 
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We also use the following inverse assumption: there exists a constant C , independent 
of h and s h , such that 


(4.83) ||s h ||,, r <C'/i , -«||s h || ?i r V s h £ P h , -1/2 < q < s < 1/2 . 

See [2] or [8] for details concerning (4.82) and (4.83). See also [15] for (4.82). 

Now, let Q h = V ft | r , he., Q h consists of the restriction, to the boundary T, of 
functions belonging to V ft . Again, for all choices of conforming finite element spaces 
we then have that C H 1 (T). We can show the approximation property: there 
exist an integer k and a constant C , independent of h and k, such that for 1 < m < k, 
0 < s < 1 and k E H m+1 (fi)| r , 


(4.84) 


inf ||k - k ft || S) r < inf ||v|| m+1 

lt h £Q h vEH m + 1 (f2),v|r=k 


This property follows from (4.79), once one notes that the same type of polynomials 
are used in as are used in V ft . We set G h = Q ft nHi(r). 

Once the approximating subspaces have been chosen we seek £ V h , p h £ S£, 
t h e g h e Q v h e v ft , e S h^ T h e and e ]R suc h that 

(4.85) a(u ft ,v ft ) + Ac(u ft ,u ft ,v ft ) + &(v ft ,/)-(v ft ,t ft ) r = A(f,v ft ) V v ft £ V h , 


(4.86) 


b(u h ,q h ) = 0 Vq h £S% 


(4.87) 


« s ft ) r - A(g ft , s ft ) r = A(b, s ft ) r V £ P h 


(4.88) (V s g ft ,V s k ft )r + (g ft ,k ft ) r + 7 ft ^k ft -ndr= (r ft ,k ft ) r V k h £ Q ft 


(4.89) 


g ft ndT = 0, 


(4.90) 


a(w h ,i/ h ) + A c(w ft , u h , v h ) + A c(u ft , w h , v h ) + b(w h , <f> h ) - {w h , r h ) T 

= A ((u ft — u 0 ) 3 , w ft ) V w h £ V h , 


(4.91) 


b(v h ,r h ) = 0 V 


(4.92) 


(v h , y h ) = 0 V £ P h . 


Note that if (4.85)-(4.92) are satished, then necessarily g ft £ G h . Also, in the right- 
hand side of (4.90), we use a notation similar to that used in the right-hand side of 
(4.70). 


42 



The operator T h £ C(Y,X h ) is defined as the solution operator for 
a(u ft ,v ft ) + &(v\/)-(v\t ft ) r = (f,v ft ) V v ft £ V\ 
b(u\q h ) = 0 Vq h eS%, 

and 

(u ft ,s ft ) r = (b,s ft ) r V s h £ P h ; 

i.e, for each f £ Y , T h f = ij) h £ is the solution of the above system of equations. 
Since T = T*, we define ( T* ) h = T h . 

We define the operator E h : G* — >■ G ft as follows. For each T £ G*, g ft = E h T if 
and only if 

(V s g ft , V s z ft ) r + (g ft , z ft ) r + 7 h J z h ■ n dr = ( r h ,z h ) r V z ft £ Q h 


and 


I g h ndT = 0. 
r 


The existence and uniqueness of a solution (g ft , 7 ft ) £ x 1R is guaranteed by the 
Brezzi theory for mixed finite element methods (see [4] or [5]) and the inequalities 


(4.93) (V s k ft ,V s k ft )r + (k ft ,k ft ) r > G||k ft ||? >r Vk h £Q h cH 1 (r) 

and 


(4.94) 


sup 

0^k'*eQ ,t 


/l 


r k ft • n dr . 


V j h £ R . 


The solution necessarily satisfies g ft £ G h . Thus the operator E h is well defined. 

With these definitions we see that (4.85)-(4.92) can be written in the form (3.1)- 

(3-3). 

By results concerning the approximation of the Navier-Stokes equations with 
inhomogeneous boundary conditions (see [15]), we obtain 

||(T — T h )f\\x —>■ 0 


as h —>■ 0, for all / = (£, rj, k) £ Y . This is simply a restatement of (H13). 

(H14) follows trivially from (H13), the fact that T is self-adjoint, and the choice 

To verify (H15), we note that the nondiscretized version of (4.93)-(4.94) certainly 
also holds, i.e., 


and 


(V s k,V s k)r + (k,k) r > G||k||f r V k £ H 1 (r) 
7 J r k ■ n dF 


sup 
o^keHqr) 


i,r 


> G| 7 | V 7 £ 1R . 


Using the Brezzi theory for mixed finite element method (see [4] or [5]), we obtain 
that 

\\(E — E h )r\\i r —>■ 0 as h 0 , 
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which verifies (H15). 

(H16) and (H17) follow from the fact that N and T are polynomials. Here we 
also use imbedding theorems and Cauchy inequalities. 

We set Z, = L 3 / 2 (£!) x {0} x {0}. For each (v,q,s) £ X = H 1 ^) x L$(fi) x 
H-!/ 2 (r) and (w,r, k) £ X = H 1 (S1) x T 2 (fl) x H _1 / 2 (r), Sobolev imbedding theo- 
rems imply that 


[^'(wM)]* • (v,g,s) 


(u- V)v + v -(Vu) T \ 

o fez, 

0 / 


([X"(U, p, t)]* • (v, q, s)) • (w ,r, k) 


(w • V)v + v • (Vw) T \ 

o fez, 

o 


and 


l 


/ 3(«i - m 0 i) 2 wi^i 


\ 


p,t)-(v,q, s)) • (w, r, k) 


\3 (u d - u 0d ) 2 w d v d 

0 

V 0 


e z , 


where d (= 2 or 3) is the space dimension. These relations verify (H18). 

From the definition of the operator K we see that K maps H^(T) into L 3 / 2 (£!) x 
{0} x H 1 (r), i.e. , K maps G into Z. Thus (H19) is verified. 

Hence, we are now in a position to apply Theorem 3.5 to derive error estimates 
for the approximate solutions of the optimality system (4.54)-(4.56), (4.70)-(4.72) and 
(4.74). It should be noted that Lemma 3.4 implies that for almost all values of A, the 
solutions of the optimality system are regular. 

THEOREM 4.10. Assume that A is a compact interval of 1R_|_ and that there exists 
a branch {(A, u(A),p(A), t(A), g(A), J/(A), <^(A), t(A)) : A £ A} of regular solutions of 
the optimality system (4.54)-(4.56), (4.70)-(4.72), and (4.74). Assume that the finite 
element spaces X h and G h satisfy the hypotheses (4.79)-(4.84). Then, there exists a 
8 > 0 and an ho > 0 such that for h < ho, the discrete optimality system (4.85)-(4.92) 
has a unique branch of solutions {(A,u h (A),p h (A),t h (A),g h (A),i' h (A),(f) h (A),T h (A)) : 
A £ A} satisfying 


(l |u( A) - u ft (A)|| 1 + \\p(A) - /(A)||o + | |t (A) - t ft ( A) 1 1 _ 1/2>r 
+ ||g(A) - g ft (A)|| lir + HA) - v A (A)||i + ||0(A) - 0 ft (A)|| o 
+ ||r(A) - r ft (A)||_ li2i r) < (5 for all A £ A. 


Moreover, 


lim (| |u( A) - u ft (A)|| 1 + \\p(A) - /(A)||o + ||t(A) - t ft (A)||_ 1/2 , r + ||g(A) - g ft (A)|| lir 
+ HA) - ^(A)!!, + ||0(A) - 0 ft (A)|| o + ||r(A) - r ft (A)||_ li2 , r ) = 0 


uniformly in A £ A. 
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If, in addition, the solution satisfies (u(A), p(A), t(A), g(A), J/(A), />(A), t(A)) E 
H m+1 (SI) x H m {Q) x H m (fi)| r x H m+1 (fl)| r x H m+1 (fi) x H m (Q) x H m (fi)| r for 
A E A, then there exists a constant C', independent of h, such that 


( ||u(A) - u ft (A)|| 1 + ||p(A) - /(A)||o + | |t (A) - t ft (A)||_ 1/2 r + ||g(A) - g ft (A)|| lir 
+ UKA) - ( A) 1 1 x + ||0(A) - 0 ft (A)|| o + ||r(A) - r ft (A)||_ li2 , r ) 

< C/j m - 1 / 2 (||u(A)|| m+ i + ||p(A)|| m + inf 1 1 v | |m 

V veH”»(n),v|r=t 


+ inf 

vgH”»+ 1 (n),v| r =g 


|m + l + || i/ (^)||m + l + 


wGH 


inf 

*(«).' 


'|r=T 


uniformly in A E A. 

Proof: All results follow from Theorem 3.5. For the last result, we also use (3.25) and 
the estimates (see, e.g., [16] or [17]) 

< Ch m (\\u\\ m+1 + \\p\\ m + inf ||v|| m ) 

veH”»(n),v|r=t 

for u E H m+1 (£2), p E H m ( Q), and t E H m (fi)| r , 


WifT^iTT 1 ~ I)(»,4 >,t)\\ Y * = \\(T h T- l -I)( V ,f,T)\\ x 

< C/l m (||l/|| m + l + IMIm + inf ll w l|m) 

weH”»(n),w|r=T 

for v E H m+1 (S2), f E H m { Q), and r E H m (fl)| r 


and 


\\(E h E~ 1 - -f)g|| lir < Ch m - 1 / 2 


vgH 


inf 

* + 1 (n).- 


'|r=g 


|m + l 


for g E H m+1 (fl)| r . 


In these estimates, the constant C is independent of h , u, p, t, g, v, <f>, r, and A. □ 
Remark. If the control g E H m+3 / 2 (S2) |r, then the exponent of h in the error estimate 
of Theorem 4.10 can be increased from (m — 1/2) to m. □ 


5. Conclusions 

We have set up an abstract framework for the analysis and approximation of 
a class of nonlinear optimal control and optimization problems. Nonlinearities can 
occur in both the objective functional and in the constraints. Within the framework we 
have defined an abstract nonlinear optimization problem posed on infinite dimensional 
spaces, defined an approximate problem posed on finite dimensional spaces, and listed 
a number of hypotheses concerning the two problems. We then have shown that 
optimal solutions exist and that Lagrange multipliers may be used to enforce the 
constraints. We then used the Lagrange multiplier rule to derive an optimality system 
from which optimal states and controls may be deduced. We then derived existence 
results and error estimates for solutions of the approximate problem. The abstract 
framework and the results derived from that framework were then applied to three 
concrete control or optimization problems and their approximation by finite element 
methods. The first involves the von Karman plate equations of nonlinear elasticity, 
the second the Ginzburg-Landau equations of superconductivity, and the third the 
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Navier-Stokes equations for incompressible, viscous flows. It is certainly possible to 
apply the abstract results that we have derived to a variety of optimal control problems 
arising in other settings. 
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